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ABSTRACT 


This  report  deals  with  pre-  and  post-test  work  on  the  Air  Force 
Structures  Test,  Project  3.4,  "Tests  on  the  Loading  of  Truss  Systems 
Common  to  Open  Frame  Structures."  The  over-all  objective  of  the  test 
was  to  determine  the  blast  loading  on  open- framed  structures,  and  to 
compare  these  loads  with  the  results  of  steady-state  wind-tunnel  tests 
on  model  structures  where  possible.  The  items  tested  were  a section  of 
A plate  girder  bridge,  a center  section  of  a truss  bridge,  the  upper 
chord  and  lower  chord  of  this  bridge,  and  a typical  isolated  structural 
member  as  represented  by  an  I-beam.  All  items  were  at  2000  ft  (intended) 
from  ground  zero  in  both  Shots  9 and  10.  Thus,  they  were  in  the  regular 
reflection  region  in  Shot  9 and  in  the  precursor  region  in  Shot  10.  In- 
strumentation consisted  of  strain  gage  measurements  taken  on  the  sup- 
porting foundation  of  the  structures. 

Thirteen  of  eighteen  strain  gages  provided  usable  records  in  Shot  9, 
and  no  usable  records  were  obtained  in  Shot  10;  even  those  records  clas- 
sified as  usable  appeared  to  have  significant  baseline  shifts.  For  the 
most  part,  the  strain  data  did  not  yield  any  interpretable  results,  and 
therefore  none  of  tne  stated  test  objectives  has  been  achieved.  In  view 
of  the  limited  and  uncertain  nature  of  the  strain  data  obtained,  it  is 
very  doubtful  if  the  test  objectives  could  have  been  realized  in  any 
event.  However,  the  anticipated  methods  of  analysis  were  unsuccessful 
for  the  data  obtained,  and  it  is  clear  that  additional  work  remains  to 
be  done  with  regard  to  the  data-reduction  problems  associated  with  net 
force  measurement  systems  such  as  employed  on  the  present  test. 

The  major  conclusion  of  the  program  resulted  from  consideration  of 
the  response  of  the  truss  bridge  section  which  sustained  a slight  perma- 
nent set  on  Shot  9 and  was  overturned  on  Shot  10.  It  was  found  that  a 
simplified  dynamic  response  analysis,  incorporating  the  (pretest)  pre- 
dicted loading  in  the  regular  reflection  region  for  Shot  9 and  a tenta- 
tive load  prediction  scheme  applicable  to  the  precursor  region  for  Shot 
10,  provided  an  adequate  estimate  of  the  damage  sustained  by  the  bridge 
in  these  two  shots.  Thus,  confidence  in  the  utilization  of  existing 
methods  for  damage  prediction  estimates  of  open-framed  structures  is 
certainly  increased  as  a result  of  this  test,  even  though  a precise 
determination  of  the  blast  loading  was  not  achieved. 
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FOREWORD 


" 


This  report  is  one  of  the  reports  presenting  the  results  of  the 
78  projects  participating  in  the  Military  Effects  Tests  Program  of 
Operation  UPSHOT-KNOTHOLE,  which  included  11  test  detonations.  For 
readers  interested  in  other  pertinent  test  information,  reference  is 
made  to  WT-782,  Summary  Report  of  the  Technical  Director,  Military  Ef- 
fects Program.  This  summary  report  includes  the  following  information 
of  possible  general  interest. 

a.  An  over-all  description  of  each  detonation,  including 
yield,  height  of  burst,  ground  zero  location,  time  of 
detonation,  ambient  atmospheric  conditions  at  detona- 
tion, etc.,  for  the  11  shots. 

b.  Compilation  and  correlation  of  all  project  results  on 
the  basic  measurements  of  blast  and  shock,  thermal 
radiation,  and  nuclear  radiation. 

c.  Compilation  and  correlation  of  the  various  project  re- 
sults on  weapons  effects. 

d.  A summary  of  each  project,  including  objectives  and 
results . 

e.  A complete  listing  of  all  reports  covering  the  Mili- 
tary Effects  Tests  Program. 


PREFACE 


In  a letter  dated  12  March  1952,  the  Air  Materiel  Command  was 
requested  by  Air  Research  and  Development  Command  to  submit  for  test- 
ing in  Operation  UPSHOT -KNOTHOLE  existing  requirements  for  a structures 
program  which  would  be  based  on  the  needs  of  the  Air  Force  for  target 
analysis  and  indirect  bomb  damage  assessment  information.  Within  the 
Air  Materiel  Command,  the  responsibility  for  designing  and  executing 
such  a program  was  delegated  to  the  Special  Studies  Office,  Engineer- 
ing Branch  of  the  Installations  Division.  The  requirements  which  were 
submitted  and  approved  became  part  of  Program  3 of  the  operation  and 
were  designated  as  Projects  3*1,  3*3,  3*^,  3*5,  3*6  and  3*26.1.  Mr. 

B.  J.  O’Brien  of  the  Special  Studies  Office  was  appointed  Project 
Officer  and,  as  such,  coordinated  and  successfully  directed  the  plan- 
ning and  operational  phases  of  five  of  the  six  projects.  Due  to  the 
similarity  in  test  objectives  involving  railroad  equipment,  the  pro- 
jects proposed  by  the  Transportation  Corps,  U.  S.  Army,  and.  the  U.  S. 
Air  Force  were  combined  into  Project  3*6  with  Lt.  Colonel  Donald  G. 

Dow,  TC,  USA,  as  Project  Officer  and  Mr.  O’Brien  as  Assistant  Project 
Officer. 

Armour  Research  Foundation  (ARF)  of  the  Illinois  Institute  of 
Technology  was  awarded  a contract  to  assist  the  Special  Studies  Office 
in  planning  and  designing  the  experiments,  and  in  analysis  and  report- 
ing of  test  results.  During  the  period  of  planning,  close  liaison 
was  maintained  with  other  interested  Air  Force  agencies,  particularly 
the  Physical  Vulnerability  Division,  Directorate  of  Intelligence, 
Headquarters,  USAF.  Many  valuable  suggestions  were  contributed  by 
Colonel  John  We It man,  USAF,  Lt.  Colonel  John  Ault,  USAF,  Messrs. 

R.  G.  Grassy  and  S.  White,  Dr.  F.  Genevese  and  others  of  that  Division, 
and  by  Mr.  Louis  A.  Nees,  Chief  of  the  Engineering  Branch,  Installa- 
tions Division,  AMC. 

Personne1  of  the  Special  Studies  Section  who  were  intimately 
connected  with  the  program  were  Mr.  Eric  H.  Wang,  Chief,  Special 
Studies  Office,  who  was  the  technical  and  scientific  monitor  for  the 
Air  Force  Program,  Mr.  Arthur  Stansel,  and  Mrs.  Maisie  G.  Ridgeway, 
secretary  to  Mr.  Wang.  Other  members  of  the  Office  who  were  associated 
with  the  program  were  Messrs.  R.  R.  Birukoff,  P.  A.  Cooley,  J.  C.  Noble, 
V and  Lts.  T.  M.  Murray  and  G.  A.  Rockwell,  USAF. 
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Most  of  the  Introduction  section  of  this  report  was  taken  from 
the  preface  of  the  Preliminary  Report,  Operation  UPSHOT -KNOTHOLE, 
Project  (3.1)  authored  by  Eric  H.  Wang  and  Bernard  J.  O’Brien. 

The  responsibility  within  the  Air  Force  for  execution  of  the  six 
projects  was  transferred  from  the  Special  Studies  Office,  Installations 
Division,  Air  Materiel  Command  to  Blast  Effects  Research,  Mechanics 
Branch,  Aeronautical  Research  Laboratory,  Wright  Air  Development 
Center,  on  15  Noveniber  195^  • 
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CHAPTER  1 


INTRODUCTION 


1.1  PURPOSE  OF  AIR  FORCE  TEST  PROGRAMS 


The  series  of  tests  conducted  by  the  Air  Force  in  Operation  UPSHOT- 
KNOTHOLE  is  part  of  a continuing  Air  Force  program  designated  as  "Deter- 
mination of  Blast  Effects  on  Buildings  and  Structures."  The  United  States 
Air  Force  is  mainly  interested  in  the  offensive  aspects  of  such  research. 

The  UPSHOT-KNOTHOLE  projects  sponsored  by  the  Air  Force  and  their 
specific  objectives  cannot  be  fully  understood  without  some  knowledge  of 
the  general  objectives  of  the  over-all  program.  The  research  results 
emanating  from  these  studies  and  experiments  conducted  by  the  Air  Force 
are  used  by  a number  of  government  agencies  to  improve  their  own  systems 
of  determining  blast  effects,  or  to  further  their  own  research. 

One  of  these  agencies  is  the  Directorate  of  Intelligence,  Head- 
quarters, USAF,  which  feeds  results  as  they  are  obtained  into  its  own 
system  of  vulnerability  classes,  thereby  making  it  possible  to  analyze 
prospective  enemy  targets  with  greater  accuracy,  and  to  recommend  the  de- 
sired ground  zero.  Another  principal  user  of  the  research  results  is 
the  Strategic  Air  Command,  which  applies  them  toward  improvement  of  an 
existing  indirect  bomb  damage  assessment  system.  The  purpose  of  this 
system  is  to  make  it  possible  to  dispense  with  the  usual  reconnaissance 
after  a strike,  using  instead  information  on  the  actual  ground  zero, 
height  of  burst,  and  yield  of  the  weapon  which  is  brought  back  to  the 
operational  base  by  the  strike  aircraft  to  determine  the  damage  inflicted. 

The  task  of  determining  the  effect  of  blast  on  various  types  of 
building  structures  and  tactical  equipment  is  a rather  formidable  one. 
However,  its  difficulty  is  somewhat  relieved  by  the  fact  that,  for  the 
offensive  purposes  in  which  the  Air  Force  is  interested,  it  is  not 
necessary  to  determine  the  effect  of  transient  loads  on  these  items  with 
the  same  accuracy  as  would  normally  be  employed  for  static  design  pur- 
poses. In  fact,  even  if  it  were  possible  to  solve  the  dynamic  problems 
satisfactorily,  Intelligence  information  would  be  far  too  sketchy  to 
furnish  the  information  necessary  to  Justify  the  use  of  an  accurate 
analysis  for  items  located  in  prospective  enemy  countries.  From  the 
experience  that  is  so  far  available  it  is  expected  that  it  will  be 
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possible  within  the  foreseeable  future  to  determine  blast  damage  within 
broad  limits  with  sufficient  accuracy  for  planning  as  well  as  for  oper- 
ational purposes. 

In  view  of  the  complex  phenomena  attending  shock  waves  emanating 
from  various  types  of  atomic  blasts  and  the  uncertainties  inherent  in 
determining  s ignif icant  parameters,  an  investigator's  first  idea  would 
be  to  obtain  solutions  through  a long  series  of  very  elaborate  and  prop- 
erly designed  full-scale  tests.  However,  neither  funds  nor  time  will 
allow  such  an  approach.  It  has  therefore  been  the  objective  of  the 
agencies  involved  to  obtain  sufficiently  accurate  results  by  judicious 
use  of  theoretical  analyses,  laboratory  tests,  high-explosive  field 
tests,  and  a small  number  of  full-scale  atomic  tests. 

Three  of  these  research  projects  have  involved  full-scale  atomic 
testing.  The  first  was  GREENHOUSE,  the  second  was  JANGLE  (the  first 
underground  hurst  of  an  atomic  weapon  to  which  an  Air  Force  structures 
program  was  subjected)  and  the  third  the  present  UKHOT-KNCHHOLE 
program. 

From  previous  analysis,  laboratory  tests,  and  full-scale  tests 
(the  latter  especially  as  conducted  in  GREENHOUSE),  methods  of  damage 
prediction  have  been  developed  by  Armour  Research  Foundation  (ARF)  and 
others.  These  prediction  methods  have  attempted  to  describe  the  char- 
acter of  the  blast  loads  acting  on  a variety  of  items.  Response  compu- 
tations based  on  the  predicted  loadings  permit,  in  turn,  an  estimate  of 
physical  damage.  However,  the  relation  between  the  deflection  or  move- 
ment of  a body  and  significant  military  damage  has  never  been  clearly 
established  except  for  extreme  cases,  e.g.,  total  destruction  or  no 
destruction.  Another  aim  of  these  tests  is,  therefore,  to  establish 
the  relationship  between  deflection  and  functional  damage.  A full-scale 
test  also  affords  an  excellent  opportunity  to  determine  scaling  check 
points  for  laboratory  tests. 

In  addition  to  the  scientific  aspects  of  the  tests,  most  of  the 
results  of  the  Air  Force  projects  can  be  used  by  other  government  agen- 
cies such  as  the  Directorate  of  Intelligence  to  furnish  "rough-and-ready" 
experimental  answers  to  the  behavior  of  various  kinds  of  structures 
under  blast.  In  many  cases  there  is  a statistically  significant  number 
of  items  involved  which,  added  to  previous  experimental  data  such  as 
those  gathered  at  Hiroshima  and  Nagasaki,  will  help  round  out  the  pres- 
ent vulnerability  picture.  In  other  cases,  mathematical  analysis  may 
have  to  rely  on  ad  hoc  information  to  furnish  parameters  which  cannot 
be  obtained  in  any  other  way. 

The  foregoing  remarks  are  designed  to  furnish  the  background  neces- 
sary for  a full  understanding  of  the  objectives  of  this  and  other  of  the 
Air  Force  projects.  The  full  significance  and  value  of  the  results  of 
each  test  will  be  realized  only  when  they  are  correlated  with  results 
of  past,  current,  and  future  analyses,  laboratory  tests,  high- explosive 
field  tests,  and  full-scale  atomic  investigations. 

1.2  SPECIFIC  OBJECTIVES 

Many  existing  items  of  military  importance  are  referred  to  as  "drag- 
type  structures."  This  nomenclature  covers  those  structures  having  a 
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relatively  small  area  as  compared  to  the  total  volume  they  encompass . 
Some,  like  truss  bridges,  are  constructed  of  small-dimensioned  com- 
ponents. Others,  such  as  buildings  whose  walls  and  roofs  axe  quickly 
destroyed  by  the  blast,  consist  of  a structural  frame  of  relatively 
small  area.  Evidence  exists  to  show  that  for  such  structures  most  of 
the  blast  damage  is  due  to  forces  acting  after  the  initial  diffraction 
period  of  the  blast  loading. 

Present  load  prediction  methods  on  drag-type  structures  are  adap- 
tations of  steady- state  concepts;  that  is,  the  net  loads  are  taken  to 
be  proportional  to  the  wind  pressures  behind  the  shock  multiplied  by 
the  drag  coefficient  of  the  structural  shape  as  determined  from  steady- 
state  measurements.  It  is  not  known  whether  the  laws  developed  for 
steady- state  drag  are  applicable  in  this  case,  or  whether  the  drag  coef- 
ficients for  various  shapes  determined  under  steady  flow  conditions  ap- 
ply. It  has  not  yet  been  possible  to  produce  transient  flow  in  wind 
tunnels  which  adequately  reproduce  field  conditions,  and  until  either 
this  or  extended  large  shock  tube  experiments  in  this  field  are  possible, 
the  collection  of  full-scale  check  data  becomes  almost  mandatory. 

The  specific  objectives  of  this  test  were: 

1.  To  determine  the  differences  between  drag  loads  due  to 
an  atomic  blast  and  those  of  an  apparently  similar  nature 
produced  in  wind  tunnels  under  steady- state  flow  con- 
ditions. 

2.  To  determine  the  ratio  of  diffraction  Jjapulse  to  drag 
impulse  in  drag- type  structures. 

3.  To  find  the  relationship  between  drag  coefficients  for 
the  same  body  when  subjected  to  transient  as  compared 
to  steady-state  flow  conditions. 

4.  To  find  the  effects  on  drag  loading  due  to  shielding 
of  component  parts. 

5.  To  determine  the  relationship  between  the  loadings  on 
drag- type  structures  which  are  scale  models. 

The  over-all  objective  of  this  test  was  to  obtain  sufficient  in- 
formation as  to  the  loading  on  drag-type  structures  so  that  damage 
estimates  can  be  made  by  means  of  analytical  techniques. 

1.3  RESPONSIBILITIES 

Armour  Research  Foundation  (ARF)  was  retained  by  the  Air  Materiel 
Command  (AMD)  of  the  United  States  Air  Force  to  carry  out  the  following 
specific  objectives  of  the  program: 

1.  Consultation  on  the  selection  of  the  test  items. 

2.  Design  of  the  test  items. 

3.  Specification  of  instrumentation  requirements. 
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4.  Location  of  the  structures  at  the  test  site. 

5.  Supervision  of  construction  of  the  test  items. 

6.  Theoretical  analyses  concerning  pretest 
predictions  of  blast  loading  and  response  of  the  test 
items  where  required. 

7.  Analysis  of  the  test  results. 

8.  Submission  of  reports  accounting  for  the  ARF 
activities  pursuant  to  the  objectives  of  the  program. 

Detailed  statements  of  the  duties  and  obligations  of  the  contracting 
parties  can  be  found  under  the  "Statement  of  Work"  in  Air  Force  Contract 
AF33(038) -30029. 

Preparation  of  the  construction  drawings  for  most  of  the  test  items 
was  subcontracted  by  ARF  to  the  firm  of  Holabird  and  Root  and  Burgee. 

A member  of  this  organization  supervised  the  actual  construction  under 
the  general  direction  of  ARF.  As-built  drawings  of  all  the  test  items 
were  prepared  by  the  Silas  Mason  Company. 

All  instrumentation  was  installed  and' operated  by  the  Ballistic 
Research  Laboratories  (BRL)  under  UPSHOT-KNOTHOLE  Project  3*28.1 
( "Structures  Instrumentation,"  WT-738) •i' 


1/  Complete  reference  concerning  publications  mentioned  in  parentheses 
may  be  found  in  the  Bibliography  at  the  end  of  this  report. 
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CHAPTER  2 

GENERAL  DESCRIPTION  OF  TEST 


2.1  TEST  ITEMS 

A series  of  five  open-framed  structures  was  included  in  both  Shots 
9 and  10.  The  structures  were  located  in  pressure  regions  where  large 
deformations  were  not  anticipated.  The  basic  open-framed  structure  was 
a duplicate  of  the  center  section  of  a through- type, open- deck,  single- 
track^  truss  bridge  (designated  as  test  item  3*4a).  Duplicates  of  the 
top  chord  assembly  ( 3.4b) , the  bottom  chord  assembly  (3.4c),  and  a single 
I-beam  from  the  latter  section  (3.4e)  were  tested.  The  fifth  test  item 
was  a section  of  an  Atchison,  Topeka  and  Santa  Pe  Railroad  Class  D plate 
girder,  through- type,  open-deck  bridge  (3.4f). 

' Each  item  was  mounted  on  specially  designed  supports,  or  sensor 

bars,  which  measured  the  reactions  induced  by  the  blast  loading.  The 
sensor  bars  were  welded  to  the  test  structures  at  one  end  and  bolted 
to  reinforced  concrete  piers  at  the  other  end.  The  design  of  the  piers 
was  such  as  to  prevent  any  gross  motion  of  the  test  structures.  Guy 
cables  anchored  to  dead-men  were  fastened  at  each  of  the  four  corners, 
just  above  the  sensor  bars,  on  the  four  larger  assemblies  tested.  Each 
cable  was  pre-tensioned  to  about  5000  lb,  and  the  ends  were  secured  with 
cable  clamps.  In  addition,  cross-bracing  cables  were  added  to  each  end 
of  the  3.4a  truss  bridge  section  in  lieu  of  end  portals.  Preshot  photo- 
graphs of  the  test  items  are  shown  in  Figs.  2.1  through  2.5.  As-built 
construction  drawings  are  available  upon  request  from  AFSWP.i/ 

2.2  INSTRUMENTATION 
2.2.1  General 

The  instrumentation  was  designed  to  measure  the  net  blast  forces 
acting  on  the  test  items.  In  view  of  the  complex  geometry  involved 
(i.e.,  many  members  of  relatively  small  dimension),  it  was  impractical 
to  obtain  direct  pressure  measurements  from  which  the  desired  net  forces 
could  be  determined.  The  system  employed  consisted  of  mounting  the  items 
on  specially  designed  supports  (see  Fig.  2.6).  Time-dependent  reactions 
17  Write  to:  Field  Command,  Armed  Forces  Special  Weapons  Project, 

Sandia  Base,  Albuquerque , New  Mexico,  Attention:  WETD. 
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transmitted  by  the  structures  to  the  supports,  or  sensors,  were  measured 
by  means  of  resistance  strain  gages.  The  strain  gages  were  mounted  on 
the  front  and  back  of  each  sensor  support  and  connected  in  series  elec- 
trically to  cancel  bending  strains.  Corresponding  sensor  bars  on  op- 
posite sides  of  the  structures  were  averaged  so  that  a total  of  four 
instrumentation  channels  per  structure  was  used  for  strain  measurements, 
with  the  exception  of  the  beam  ( 3.4e)  which  used  only  two  channels. 

The  gages  were  numbered  consecutively  in  the  downstream  direction,  i.e., 
strain  gage  SI  was  on  the  upstream  sensor  and  strain  gage  34  was  on  the 
furthest  downstream  sensor. 

The  strain  data  represent  the  output  of  a so-called  net  force 
measurement  system.  It  should  be  recognized,  however,  that  the  strain 
data  as  such  do  not  give  a direct  measure  of  the  applied  blast  loading, 
but  rather  the  response  of  the  structure-sensor  system  to  this  loading. 
In  effect,  the  test  structures  are  being  utilized  as  dynamometers,  and 
a suitable  analysis  of  the  strain  data  is  required  in  order  to  determine 
the  desired  information. 

2.2.2  Strain  Measurements 

Standard  SR4  strain  gages  were  used  in  a four  active  arm  bridge 
configuration  to  measure  axial  strain.  The  output  of  each  bridge  was 
fed  into  a Webster-Chicago  recording  system  through  a coupling  unix  and 
recorded  on  magnetic  tape. 

The  calibration  of  the  strain  gages  was  accomplished  electrically 
by  shunting  the  proper  arm  of  each  gage  installation  with  an  accurately 
known  resistance  to  simulate  actual  strain. 

Complete  details  of  the  strain  gage  installations  are  contained 
in  the  final  report  of  UPSHOT-KNOTHOLE  Project  3.28.1  ("Structures  In- 
strumentation," WT-738) • 

2.2.3  Instrument  Records 

The  Ballistic  Research  Laboratories  (BRL)  handled  all  of  the  in- 
strumentation oh  the  program.  The  output  of  the  strain  gages  was  re- 
corded initially  on  magnetic  tape  and  later  played  back  onto  oscil- 
lographic paper.  The  records  in  this  form  exhibit  characteristics  which 
made  them  undesirable  for  purposes  of  interpretation  and  comparison 
(e.g.,  the  ordinate  scale  is  non-linear).  For  that  reason,  all  of  the 
records  were  converted  into  linear  form. 

BRL  reduced,  calibrated,  and  plotted  to  linear  scales  all  of  the 
strain  records.  The  records  as  presented  to  ARF  were  in  the  form  of 
plots  made  up  of  the  points  at  which  the  records  were  read.  ARF  was 
responsible  for  fairing  curves  through  these  points.  BRL  also  submitted 
tabulated  listings  of  the  points,  as  well  as  copies  of  the  original 
playbacks . 

The  numerical  data-reduction  scheme  employed  on  the  strain  re- 
cords (see  Section  5. 1.1)  required  that  the  data  be  given  in  equal  time 
intervals.  Accordingly,  BRL  made  a separate  linearization  of  selected 
strain  records  at  equal  time  intervals  of  about  4 ms,  and  supplied  these 
results  on  IBM  punched  cards. 
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Actual  G.Z.,  Shot  10;  H • 524' 


Actual  6.Z.,  Shot  9;  H • 2423' 


Not*;  All  item*  included  in  both  Shot*  9 and  10. 
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Fig.  2.7  Location  of  Structures  at  Test  Site 


2.3  LOCATION  OF  TEST  STRUCTURES 

The  location  of  the  structures  at  the  test  site  is  shown  in  Fig. 
2.7.  A summary  of  the  field  conditions  including  ground  range  from 
intended  and  actual  ground  zero  and  measured  overpressure  levels  is 
given  in  Table  4.1.  The  maximum  misorientation  resulting  from  the  bomb- 
ing error  was  about  22  degrees  in  Shot  9 and  about  5 degrees  in  Shot  10. 
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CHAPTER  3 


PRETEST  CONSIDERATIONS 


3.1  SELECTION  OF  TEST  ITEMS 

The  pretest  work  on  this  program  is  contained  in  the  Planning  Pro- 
gram for  Air  Force  Structures  Tests,  Final  Pretest  Report  (Part  VIII, 
"Tests  on  Open-Grid  Structures"). 

After  consultation  with  the  Air  Force,  the  full-scale  bridge  sec- 
tions (3.4a  and  3.4f)  were  selected.  The  specific  designs  were  chosen 
primarily  because  construction  drawings  were  readily  available. 

The  number  and  disposition  of  the  test  structures  was  chosen  in 
order  to  achieve  the  stated  objectives  with  the  fewest  items  possible. 

It  was  felt  that  the  objectives  dealing  with  transient  versus  steady- 
state  drag  effects  could  be  attained  from  the  loading  data  on  all  five 
of  the  test  items.  The  effects  of  shielding  should  be  apparent  from 
a comparison  of  the  drag  forces  on  the  component  sections  with  the 
forces  on  the  entire  bridge  section. 

It  was  planned  to  compare  the  test  data  with  the  results  of  lab- 
oratory tests  whenever  possible.  With  respect  to  shielding  effects  it 
was  hoped  that  the  test  would  serve  to  either  confirm  or  suggest  modi- 
fications of  standard  engineering  approaches  to  this  problem.  For 
example,  it  is  common  practice  to  assume  that  the  drag  forces  due  to 
winds  acting  on  the  front  or  upstream  truss  of  bridges  axe  twice  those 
acting  on  the  rear  truss.  Steady-state  wind-tunnel  tests  conducted  at 
Massachusetts  Institute  of  Technology  ("Criteria  for  Collapse  Damage 
To  Structures  in  Relation  to  Weapon  Yield,  Height  of  Burst  and  Distance") 
have  dealt  with  this  type  of  shielding  as  a function  of  the  distance 
between  trusses. 

3.2  LOAD  PREDICTIONS 


The  Planning  Program  Final  Pretest  Report  (Part  VIII,  "Tests  on 
Open-Grid  Structures")  contains  a detailed  study  of  the  anticipated 
blast  loading  on  truss  systems  situated  in  the  Mach  and  regular  reflec- 
tion regions.  This  work  is  based  on  steady-state  drag  concepts  and 
does  not  account  for  shielding  effects  of  multiple  truss  systems.  A 
summary  of  the  pertinent  results  is  presented  in  this  section. 
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3.2.1  Loading  In  the  Mach  Region 


A number  of  relationships  must  be  considered  with  reference  to 
shock  phenomena  in  the  Mach  reflection  region.  The  shock  strength, £ 
is  defined  as  the  ratio  of  absolute  pressures  across  the  Incident  shock 
front,  and  is  given  by 


c - 1 + (pa  /P0) 


(3.1) 


where  pa  is  the  initial  overpressure  (i.e.,  pressure  in  excess  of 
atmospheric)  and  PQ  is  the  ambient  atmospheric  pressure. 

The  overpressure-time  variation  of  the  blast  wave  is  approximated 
by  the  relation, 

-ctAo  . ..  . . 


pcr  ^ “ p«y  (°)e 


(1  - t/tj 


(3.2) 


where  to  is  the  duration  of  the  positive  phase  of  the  blast  wave,  and 
the  empirical  factor  c is  usually  taken  as  unity.!/  The  nominal  drag 


pressure,  p,(t) , is  given  approximately  by. 


pd(t)  » Pd(0)e 


-2ct/tQ 


(i  - t/t 0y 


(3.3) 


where 


Pd(°) 


2.5PO-  (0) 

7*0  + pcr  (0) 


(3.4) 


The  shock  front  velocity,  U,  is  given  by 


U = “T  \|  1 + 6 

NT 


= 422  \|  1 + 6(  , ft/sec  (3.5) 


where  cq  Is  atmospheric  sound  velocity,  evaluated  above  for  standard 
atmospheric  conditions. 

When  a plane  shock  front  reflects  from  a rigid  plane  surface, 
the  pressure  behind  the  reflected  wave,  pr,  is  given  by 


Pr  * P<T(0) 


8£  + 6 

■frr 


(3.6) 


The  factor  c is  known  to  depend  on  scaled  height-of -burst  and 
ground  range.  For  test  conditions,  c = 1.3  represents  a slightly 
better  approximation  to  the  free-stream  pressure  wave  form  for  both 
Shots  9 and  10  than  does  c = 1.  Plots  of  c as  a function  of 
height-of -burst  and  ground  range  are  included  in  the  final  report 
on  UPSHOT-KNOTHOLE  Project  3.1  ("Tests  on  the  Loading  of  Building 
and  Equipment  Shapes,"  WT-721) . 
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The  loading  scheme  for  open-grid  systems  i6  an  outgrowth  of 
methods  originally  presented  in  the  Operation  GREENHOUSE  Final  Report, 
Project  3.3  which,  in  turn,  were  developed  from  shock  tube  and  wind- 
tunnel  data  (Princeton  University  Technical  Reports  and  wind-tunnel 
studies  by  Chien,  Feng,  Wang,  and  Siao  and  Eiffel)  on  simple  geometric 
6hapes.  The  loading  on  these  structures  consists  of  two  phases:  the 
initial  diffraction  phase* which  depends  primarily  on  characteristic 
dimensions  of  the  structure,  and  the  subsequent  drag  phase  which  depends 
on  the  geometry  of  the  structure  as  well  as  on  the  blast  wave. 

For  a single  truss  element  (e.g.,  an_i6olated  beam)  the  c'  ■■ac- 
teristic  dimension  or  clearing  distance,  h,  refers  to  the  half -width 
of  the  member.  For  a composite  of  such  members,  the  situation  is  com- 
plicated. It  is  believed  that  an  equivalent  clearing  distance  for  the 
entire  assembly,  h,  can  be  taken  as  a weighted  average  of  the  clearing 
times  for  individual  members  (see  Operation  GREENHOUSE,  final  report). 


that  is. 


(3.7) 


where  hi  is  the  clearing  distance  for  the  ith  member  and  Ai  is  the 
frontal  area  of  that  member.  The  denominator  of  Equation  3*7  represents 
the  total  frontal  area  of  the  truss. 

The  force  per  unit  frontal  area  of  the  first  (i.e.,  upstream)  truss 
in  a structure  similar  to  3*1*®  is  shown  symbolically  in  Fig.  3»1»  The 
diffraction  loading  pattern  is  modeled  after  the  loading  on  the  two-di- 
mensional Princeton  1:2  block  (Air  Force  Planning  Program,  Final  Pretest 
Report,  Part  VIII,  "Tests  on  Open-Grid  Structures").  The  average  pres- 
sure on  the  front  surfaces  rises  instantaneously  to  the  reflected  pres- 
sure, pr,  and  then  drops  to  side-on  plus  drag  pressure, 
p (t)  + CdfPd(t),  in  an  additional  7h/U  time  units.  A combined  drag 
coefficient,  C<^  of  2 i6  selected.  This  leads  to  a value  of  the  drag 
coefficient  of  the  front  surface  of  Caf  = 1,  which  is  the  maximum  value 
that  can  be  obtained  under  steady-state  conditions.  The  drag  coefficient 
for  the  back  surface  is  taken  to  be  C,^  = -1.  Therefore,  the  net  drag 
force  on  the  truss  is  2pa(t). 

The  loading  pattern  on  the  second  (i.e.,  downstream)  truss  is 
similar  to  the  above,  with  minor  exceptions  introduced  by  the  choice  of 
drag  coefficients.  The  value  of  C<}  is  taken  to  be  1 (i.e.,  half  that 
of  the  first  truss)  in  order  to  accomt  for  shielding  effects.  This 
choice,  although  arbitrary,  is  not  in  contradiction  with  experimental 
data  ("Wind-Tunnel  Studies  of  Pressure  Distribution  on  Elementary  Build- 
ing Forms,"  State  University  of  Iowa).  To  arrive  at  this  value  of  C<j., 
the  front  and  rear  surface  coefficients  are  chosen  to  be  = 0 and 

= “I  • 

On  order  to  compute  the  net  loading  on  the  combined  trusses,  it 
is  necessary  to  introduce  a time  delay  of  L/u  in  the  loading  of  the 
second  truss  (where  L is  the  distance  between  trusses).  This  is  the 
length  of  time  relative  to  the  first  truss  required  for  the  shock  to 


reach  the  second  truss.  Figure  3.2  shows  the  loading  (i.e.,  force  per 
unit  frontal  area)  on  the  second  truss.  The  net  loading  on  the  struc- 
ture can  be  found  by  superposing  (graphically)  the  component  loadings 
on  each  truss. 

3.2.2  Loading  in  the  Regular  Reflection  Region 

The  general  methods  for  computing  loading  on  solid  structures  in 
the  regular  reflection  region  has  been  described  in  the  Final  Pretest 
Report,  Part  V ("Regular  Reflection  on  Cubical  Structures").  Since 
open-grid  structures  in  general  do  not  obstruct  the  shock  wave  with 
solid  boundaries  of  large  cross-sectional  area,  and  since  there  are  no 
re-entrant  corners  in  this  type  of  structure,  the  diffraction  phase 
will  last  for  only  a comparatively  short  time  and  the  average  loading 
will  not  build-up  to  as  high  a value  as  in  the  Mach  region.  This  is 
one  of  the  prime  differences  in  the  loading  of  solid  and  open-grid 
structures  in  the  regular  reflection  region. 

There  exist  certain  similarities  and  differences  between 
the  loading  on  structural  components  located  in  the  regular  reflection 
and  Mach  regions.  Among  the  similarities  is  the  fact  that  the  loading 
is  probably  two-dimensional  in  character,  and  hence,  the  concept  of  an 
average  h,  (Equation  3*7)  can  be  retained.  The  drag  coefficients  are 
also  assumed  to  be  the  same  for  both  cases.  In  addition,  the  ground  re- 
flected pressure,  pre(t),  is  taken  to  have  the  same  time-dependence  in 
the  regular  reflection  region  as  does  the  side-on  pressure,  p^.  (t),  in 
the  Mach  region. 

Among  the  major  differences  in  the  loading  is  the  manner  of 
reaching  pseudo- steady- state  (drag  phase)  in  regular  reflection.  Here 
the  structure  is  struck  by  two  separate  shocks,  the  incident  and  ground 
reflected  shock,  and  the  load  builds  up  in  two  distinct  time  intervals. 
This  phase  delay  is  a function  of  the  height  of  the  structure  above 
ground  a,  the  angle  of  incidence  of  the  blast,  Q , and  the  shock  front 
velocity  U.  In  addition  to  this  difference,  the  shock  waves  do  not 
envelope  the  structure  immediately,  but  rather  an  additional  phase  de- 
lay is  encountered  due  to  the  inclination  of  the  shock  fronts.  This 
inclination  results  in  a finite  time  for  the  shock  wave  to  sweep  the 
total  height,  H,  of  the  structure. 

In  considering  the  diffraction  loading  on  the  front  surface,  it 
is  convenient  to  consider  first  the  limiting  cases  of  head-on  reflection 
( Q = ir/2)  and  glancing  incidence  ( a = 0).  The  peak  pressure  re- 
sulting from  the  incident  shock  occurs  at  the  time  when  the  shock  has 
covered  the  entire  height,  H,  of  the  truss;  namely,  at  t = H cosa/U; 
this  reduces  to  zero  and  H/U,  respectively,  for  the  two  limiting  cases 
mentioned  above.  For  head-on  reflection  with  the  truss,  the  peak  aver- 
age pressure  at  this  time  is 

P - V+8i  *«■  <°>  - £ > P<r  <°>  at  ‘ ■ a-TT£  - 0 <3-6> 
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while  for  glancing  incidence  it  is 


P 3 Pg.  (0)  at,  t = 


H cos  a H 


(3.9) 


For  an  arbitrary  incident  angle,  G , the  reflected  pressure 
must  be  a function  of  8 , the  shock  strength  of  the  free-stream  inci- 
dent wave,  £ , and  the  ratio  of  the  height  of  truss,  H,  to  the  relief 
distance,  3E.  By  dividing  the  truss  into  infinitely  small  segments 
and  summing  the  average  pressure  on  each  segment  at  the  time 
t " H cos  a/lJ,  one  obtains  for  the  average  pressure  above  side-on 


Pa  (0)(C  - 1) 
2(2-77  ) 


(3-10) 


H cos  a 


P_(0)(c-1) 

-£■  ■ ^ (3.11) 


rj  m H cos  a > lm 
3b 

In  both  instances  the  factor  C is  referred  to  as  the  modified  oblique 
reflection  coefficient.  According  to  Equations  3.8  and  3*9,  the  reflec- 
tion coefficient,  C,  is  equal  to  Q (£  ) for  a = ir  / 2 and  equal  to 
unity  for  a = 0.  If  the  average  pressure  behind  the  shock  front  in 
oblique  reflection  is  plotted  against  the  angle  of  incidence,  the  re- 
sulting curve  is  approximately  a straight  line.  On  this  basis  an  em- 
pirical relationship  for  the  reflection  coefficient  as  a function  of 
the  incidence  angle  is  taken  to  be, 


■¥■  & <« > - *]  ♦ 1 


(3.12) 


When1?  £ 5,  the  average  pressure  rises  very  rapidly  to  the  value 
given  by  Equation  3.11.  This  rise  time  can  reasonably  be  neglected,  in 
which  case  the  average  pressure  at  t = 0 is  taken  to  be 


P,,.  (0)(C  - 1) 


for  V 


The  average  pressure  then  rises  linearly  from  this  value  to  p^  + p^  (0) 
at  t = H cos  a /U.  Now  for  77  < 1,  the  average  pressure  rises  linearly 
from  zero  at  t = 0 to  p.  + Pg-To)  at  t * H cos  <*  /U,  where  here  pj^ 
is  given  by  Equation  3. 10.  For  i < 77  <5,  the  average  pressure  at 
t = 0 is  obtained  by  interpolating  linearly  for  77  between  the  values  for 
T)  = 1 and  V 3 5* 
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The  pressure  build-up,  P2,  due  to  reflection  when  the  ground- 
reflected  wave  covers  the  truss,  can  be  ejqpressed  In  a similar  fashion 
where  now,  however,  pa  must  be  replaced  by  the  ground  reflected  pres- 
sure Pj.e.  The  pressure  P2  initiates  at  a phase  delay  of  2a  cos  a /U 
time  units  (a  is  the  distance  from  ground  to  the  bottom  of  the  truss, 
see  Pig.  3.3)  to  allow  for  the  reflected  shock  to  reach  the  structure. 
The  pseudo-steady-state  pressures  on  the  front  surfaces  are  taken  tobe 
Va  (0)  + ccLfPdl(°)  8ln  a for  the  Incident  wave,  and  p^O)  + cdfpdr^°) 
for  the  reflected  wave.  The  zero  times  are  again  used,  since 
the  time  interval  during  which  this  load  acts  is  very  short  and,  hence, 
the  loads  have  approximately  their  initial  value. 

Expressions  for  the  reflected  pressure,  Pre(t),  and  the  nominal 
drag  pressures,  Pai(t),  and  p^-Ct) , were  derived  in  the  pretest  re- 
port. They  are  as  follows  :£/ 

-t/t 

Pre(t)  3 Pre(°)e  (l  - t/tQ)  (3*13) 


and. 


where 


Pdr(t)  * 


pdr(0) 


1 + 


P*  <0) 


10  P, 


0 J 


-2t/t 

Pdr(0)e  0 (1  - t/tQ)2  (3.14) 


lop  2(p)  [ 6pff(0)  + 7P0] 

>„  *\  (0)]  2 


4 2 „ 
sin  a 


(3.15) 


is  the  drag  pressure  behind  the  ground-reflected  shock  front  parallel 
to  the  ground.  Also 

2.5  p<r2(0) 

- 7F  ♦,  To)  (cf-  Eq-  3-4>  <3'l6) 

O <T 

where  pdl(0)  is  nominal  drag  pressure  behind  the  incident  shock  front 
normal  to  the  incident  shock  front.  Hence,  the  component  normal  to  the 
structure  is  given  by  Pai(O)  sin  <*•  • 

The  diffraction  loading  on  the  front  of  the  second  truss  is  as- 
sumed to  be  identical  to  the  loading  on  the  front  of  the  first  truss  if 
the  proper  time  delay  for  the  shock  to  reach  the  rear  truss  is  taken 
into  account,  i.e.,  a time  delay  of  t = L sin  Ct/u.  The  drag  coefficients 
for  the  front  surfaces  are  taken  the  same  as  in  Mach  reflection,  i.e., 

» 1 on  the  front  surface  of  the  first  truss,  and  = 0 on  front 
surface  of  the  second  truss. 

2/  The  empirical  exponent  c which  determines  the  wave  form  of  the 
pressure  is  chosen  as  unity  throughout  this  section,  see  footnote 
on  page  25. 
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The  loading  on  the  downstream  surfaces  also  build  up  in  two  steps 
since  two  waves  are  involved.  The  build  up  time  for  each  wave  is  given 
by 


H cos  q 'Jh 

U U 

m 


sin  a 
U 


(3.17) 


where 


a = limiting  angle  for  regular  reflection  ( a < ctex^/ 

l.e.,  Mach  reflection  occurs  whenever  aext  is  exceeded) 
Um  =*  the  equivalent  Mach  shock  front  velocity 

= 422  ^ 1 + 6i  2 

H = depth  of  the  truss  in  the  direction  of  flow. 

The  term  i sin  a/U  is  the  time  required  for  the  wave  to  reach 
the  back  surface  of  a truss.  The  time  between  the  two  buildups  (inci- 
dent and  reflected  waves)  is  given  by  2a  cos  a/U,  which  is  the  time 
required  for  the  ground-reflected  shock  to  reach  the  back  surface. Again 
the  pseudo-steady  state  values  are  given  by  p^O)  + Cd- pdl(0)  sinOt 
for  the  incident  wave  and  Pre(°)  + cdbpdr(°)  for  the  reflected 
wave.  The  drag  coefficients  are  taken  the  same  as  for  Mach  reflection, 
namely  Cdb  = -1  for  the  back  surfaces  of  both  the  first  and  second 
trusses. 

Loading  on  front  and  rear  surfaces  of  the  first  truss  is  shown 
symbolically  in  Figs.  3.3  and  3.4  and  Tables  3.1  and  3.2.  The  loadings 
on  the  rear  truss  are  similar  to  those  on  the  front  truss.  The  time 
details  of  the  diffraction  phase  are  the  same  but  initiate  at  a time 
L sin  a /U  later  than  on  the  front  truss,  where  L is  the  distance 
between  trusses.  The  pseudo-steady-state  pressure  on  the  front  of  the 
second  truss  is  Pre(t)»  rather  than  pre(t)  + C^p^t)  , while  on 
the  rear  of  both  trusses  it  is  p-e(t)  - p^t)  . Hence,  the  total 
drag  on  the  front  truss  is  2pdr(t),  while  on  the  rear  truss  it  is 
p^Tt),  which  accounts  for  the  shielding  effect  of  the  front  truss. 

3.3  NET  FORCE  MEASUREMENTS 

3.3.I  Sensor  Design 

An  important  part  of  the  planning  phase  of  this  program  was  the 
development  of  a method  for  direct  measurement  of  the  net  blast  forces 
acting  on  the  test  structures.  Ideally,  a direct  measurement  system 
would  be  one  whose  output  is  either  proportional  to  the  net  blast  loading 
over  the  entire  structure,  or  one  whose  output  can  be  converted  to  net 
loads  in  a straightforward  fashion.  This  requirement  for  a net  force 
measurement  system  is  meant  to  exclude  pressure  gage  instrumentation, 
at  least  in  the  present  case. 


MASHED 


* 

Table  3.2  - Loading  on  Rear  Surface  of  First  Truss, 
Regular  Reflection  Region 


* For  the  back  surface  of  the  second  truss  all  times  are  increased  by 
L sin  a /U>  and  = -1. 


While  there  are  certain  types  of  structures  for  which  net  forces 
over  the  entire  structure  can  be  adequately  determined  by  means  of  com- 
bining and  averaging  individual  pressure  gage  data,  there  are  many 
structures  where  such  averaging  is  either  subject  to  considerable  un- 
certainty or  is  totally  impractical.  There  are  also  cases  where  the 
presence  of  pressure  gages  may  sufficiently  alter  the  flow  conditions 
around  the  test  items  so  as  to  rule  out  this  means  of  force  measurement. 
In  the  present  case  both  of  these  objections  can  be  raised  to  the  use 
of  pressure  gage  instrumentation.  Thus,  an  alternate  force  measurement 
system  had  to  be  devised. 

In  the  present  setup,  each  test  structure  was  utilized  in  the 
sense  of  a dynamometer.  That  is,  the  structure  and  the  sensor  arrange- 
ment as  a whole  constituted  the  gaging  device.  The  output  of  the  strain 
gage  system  (i.e.,  the  strain  in  the  sensor  bars)  was  thus  dependent 
upon  the  dynamic  characteristics  of  the  structure-sensor  system  as  well 
as  upon  the  applied  blast  force. 

In  order  for  the  strain  data  to  be  utilized  as  a direct  measure 
of  the  applied  blast  loadings,  it  is  necessary  that  the  frequency 
response  of  the  test  items  be  considerably  greater  than  that  of  the  sen- 
sor supports  alone  and,  in  turn,  that  the  fundamental  period  of  the 
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system  as  a whole  be  sufficiently  less  than  the  significant  time  details 
of  the  applied  blast  loading.  In  this  case,  the  response  is  that  of  a 
single-degree-of -freedom  system  (i.e.,  a rigid  structure  on  flexible 
supports),  and,  with  suitable  damping,  the  strain  data  could  be  inter- 
preted much  as  is  pressure  data  - simply  by  Inspection.  Unfortunately, 
it  was  not  possible  to  satisfy  these  conditions  for  any  of  structures 
tested.  The  frequency  response  of  the  sensor  bars  was  also  limited  by 
the  requirement  that  a sufficiently  great  strain  be  produced  to  yield 
a reliably  measurable  signal  from  the  strain  gages  throughout  most  of 
the  loading  period.  Also,  the  sensor  arrangement  wfs  designed  on  the 
basis  of  there  being  no  provision  for  pretest  field  calibration  of  the 
entire  setup  (and  there  was  none) . This  necessitated  a longer  sensor 
bar  than  might  otherwise  have  been  utilized,  since  it  was  imperative 
that  the  measured  surface  strain  be  an  adequate  measure  of  the  actual 
force  in  the  bar  ..2/  This,  in  turn,  further  reduced  the  frequency  re- 
sponse of  the  setup. 

Details  of  the  sensor  bar  construction  are  shown  in  Fig.  2.6, 
and  in  the  as-built  construction  drawings.  Briefly,  the  end  support 
consisted  of  a rigid  right-angle  frame  of  welded  construction  attached 
at  the  four  corners  of  the  bridge  sections, and  at  the  two  ends  of  the 
beam.  On  the  bridge  sections  the  sensor  bars  were  4-1/4  in.  square  and 
approximately  3 ft  long  between  welds  (i.e.,  the  unsupported  sensor  bar 
length);  on  the  beam  ( 3.4e)  the  bars  were  l-l/2  in.  square  and  approxi- 
mately 15  in.  long.  High-strength  steel  was  used  for  all  sensor  bars. 

The  geometry  of  the  supports  was  such  that  each  sensor  bar  could  be 
considered  as  a two-force  member.  This  would,  of  course,  be  true  if 
the  frame  were  ideally  pin-connected ,zJ  but  a rigid  frame  was  more  de- 
sirable from  a construction  viewpoint. 

3.3*2  Interpretation  of  Strain  Measurements 

The  final  sensor  design  was  such  that  the  measured  forces  in  the 
sensor  supports  would  most  likely  be  influenced  by  the  dynamic  character- 
istics of  the  entire  test  structure  and,  thus,  would  not  be  a direct 
measure  of  the  applied  blast  force.  This  situation  necessitated  a data- 
reduction  system  that  was  capable  of  taking  this  effect  into  account  and 
providing  the  desired  net  applied  forces.  In  essence  the  problem  here 
is  the  inverse  of  the  conventional  dynamic  response  problem.  That  is, 
the  usual  problem  is  "given  the  dynamic  system  and  the  input,  find  the 
output";  here  the  problem  is  "given  the  dynamic  system  and  the  output, 
find  the  input." 


17  That  is,  the  sensor  bar  had  to  be  sufficiently  long  so  that  the 

stress  distribution  at  the  cross  section  being  gaged  was  essentially 
uniform. 

4/  An  analysis  of  the  maximum  error  introduced  by  the  assumption  of  pin 
~ ends  was  carried  out  in  the  final  pretest  report ("Planning  Program 
for  Air  Force  Structures  Test, "Part  VIII,"Tests  on  Open -Grid  Struc- 
tures"). The  error  in  the  magnitude  of  the  force  was  found  to  be 
less  than  6 per  cent,  and  about  1.5  degrees  in  direction  relative 
to  the  centerline  of  the  bars. 
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The  solution  to  the  inverse  problem  is  generally  the  more  com- 
plicated but  in  principle  can  be  achieved,  at  least  for  linear  dynamic 
systems,  provided  the  response  of  the  system  to  a given  input  is  known. 
However,  the  present  problem  is  made  more  difficult  by  the  fact  that  the 
characteristics  of  the  dynamic  system  itself  are  not  known  with  any  cer- 
tainty. That  is,  the  system  is  sufficiently  complicated  that  the  response 
(i.e.,  the  force  in  the  sensor  bars)  to  a known  input  cannot  be  computed 
with  sufficient  accuracy  to  Justify  a solution  to  the  inverse  problem  by 
standard  methods.  (The  beam  (3.4e)  may  be  an  exception  to  this  state- 
ment . ) 

A numerical  data-reduction  scheme  referred  to  as  transient 
analysis"  was  developed  during  the  course  of  the  pretest  work  and  was 
believed  to  be  adequate  for  the  needs  of  this  program.  In  principle, 
the  method  of  transient  analysis  is  applicable  to  the  determination  of 
the  time  variation  of  the  input  function  (the  magnitude  of  the  input  can 
only  be  determined  to  within  an  arbitrary  scale  factor)  from  a known  out- 
put function,  even  when  the  equations  of  motion  of  the  system  are  not 
known  in  detail.  The  method  is  restricted  to  dynamic  systems  character- 
ized by  linear  second-order  differential  equations  with  constant  coef- 
ficients, and  input  functions  which  can  adequately  be  represented  by 
finite  sums  of  exponential  functions  (i.e.,  "e"  functions)  all  of  whose 
coefficients  and  exponents  are  real  numbers.  In  effect,  this  last  re- 
quirement limits  the  input  functions  to  those  having  no  harmonic  compo- 
nents, which  seemed  a reasonable  enough  assumption  in  the  case  of  blast 
loading  functions.  A further  restriction  is  that  the  forcing  function 
components  associated  with  the  various  modes  of  the  system  are  either 
zero  or  timewise  proportional  ( see  Appendix  A) . 

The  restrictive  assumptions  which  permit  the  valid  use  of  trans- 
ient analysis  as  a technique  for  data  reduction  in  the  present  applica- 
tion were  considered  to  be  satisfied  during  the  pretest  planning  period. 

As  discussed  in  Chapter  5,  the  initial  application  of  transient  analysis 
to  the  strain  data  was  largely  unsuccessful.  However,  a brief  second  ap- 
plication gave  more  encouraging  results.  The  development  and  critical 
discussion  of  transient  analysis  is  presented  in  Appendix  A of  this  re- 
port. 


3.3.3  Relationship  Between  Strain  and  Forces 

The  relationship  between  strain  and  total  force  in  each  sensor 
bar  is  given  by, 

S =*  A E « (3*18) 

8 

where 

S = force  in  sensor  bar,  lb 

€ s measured  strain  in  sensor  bar,  in. /in. 

E = Young's  modulus  for  the  sensor  bar  material,  psi 

2 

A = cross-sectional  area  of  sensor  bar,  in. 
s 
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Figure  3.5  shows  a free-body  diagram  of  a typical  test  structure 
and  sensor  assembly  applicable  to  all  items  with  the  exception  of  the 
3.4e  beam  which  was  supported  on  only  one  sensor  assembly  per  end.  With 
reference  to  Fig.  3.5,  the  incident  blast  loading  is  assumed  to  be  span- 
wise  symmetrical,  and  is  represented  by  horizontal,  F^,  and  vertical 
Fv,  time-dependent  force  components.  The  sensor  bar  forces  represent 
the  total  force  in  corresponding  sensors  on  opposite  ends  of  the  test 
structure  and,  as  such,  are  assumed  to  be  numerically  equal  to  twice 
the  measured  force  in  a single  bar  as  given  by  Equation  3.18;  hence 
the  factor  2 indicated  in  Fig.  3.5.  The  length  M is  the  distance  from 
the  vertex  of  the  sensor  assembly  to  the  center  of  pressure  of  the  test 
structure;  the  length  N is  the  distance  between  the  sensor  assemblies 
in  the  direction  of  flow. 

If  the  bottom  chord  members  are  assumed  to  be  effectively  rigid 
and  transmit  one-half  the  applied  horizontal  load  to  each  sensor  as- 
sembly, and,  further,  if  the  sensor  assembly  is  assumed  to  be  pin- 
connected  (see  Section  3.3.I)  the  determination  of  the  forces  in  the 
individual  sensor  bars  becomes  a statically  determinate  problem.  From 
a static  point  of  view,  the  horizontal  and  vertical  components  of  the 
blast  loading  are  then  given  in  terms  of  the  sensor  bar  forces  at  each 
instant  of  time  by. 


Fh  = 2 {2  (S1  - S2)  = 2/2  (S  - S4),  lb" 
Fy  = -2  \J~2  (S2  + S ) - -2/2  (S1  + S4),  lb  ^ 


(3.19) 


where  the  sensor  bar  forces,  . ..S^,  are  assumed  positive  in  tension 
as  shown  in  Fig.  3.5;  Fh  is  positive  directed  downstream  and  Fy  is 
positive  directed  downward. 

Note  that,  on  the  basis  of  this  being  a statically  determinate 
problem,  the  applied  loading  can  be  determined  from  the  measurements 
taken  at  any  one  sensor  assembly. 

The  3*^e  beam  is  supported  by  only  one  sensor  assembly  per  end, 
and  the  distance  M is  approximately  zero.  Thus,  for  this  case  the 
horizontal  and  vertical  force  components  are  given  by, 


Fh  = /2  (Sx  - S2),  lb  ' 
Fy  - [2  (S1  + S2),  lb  ^ 


(3.20) 


The  applied  forces  can  be  expressed  in  terms  of  the  measured 
strains  by  means  of  the  above  equations.  That  is,  for  structures  3.4a, 
b,  c and  f 


Fh  - 1532  ( - «2)  = 1532  ( «3  - «4),  lb  ' 
Fy  = 1532  ( <2  + «3)  = 1532  («!+%)»  lb. 


(3.21a) 
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CHAPTER  4 

EXPERIMENTAL  RESULTS 


4.1  FIELD  CONDITIONS 


All  of  the  test  structures  were  located  in  the  regular  reflection 
region  of  Shot  9,  and  within  the  precursor  region  of  Shot  10.  The  field 
conditions  are  summarized  in  Table  4.1.  As  indicated  in  the  table,  the 
bombing  error  n Shot  9 resulted  in  misorientations  up  to  22  degrees . 

4.2  PHYSICAL  IMAGE  TO  TEST  STRUCTURES 


The  blast  caused  little  damage  to  the  test  items  in  Shot  9.  The 
damage  that  did  occur  consisted  of  a small  permanent  set  of  approxi- 
mately 3 in.  at  the  top  of  the  truss  bridge  section  (3.4a)  and  the  crack- 
ing of  two  sensor  bars  at  the  '-eld  of  the  plate  girder  bridge  section 
(3.4f).  In  addition,  guy  cables  became  loose  - apparently  due  to  the 
cable  clamp  slipping.  All  cables  were  tightened  and  cracks  welded  be- 
fore Shot  10.  Two  additional  clamps  were  installed  on  each  cable,  mak- 
ing five  at  each  end.  The  permanent  deformation  of  structure  3.4a 
could  not  be  corrected  before  Shot  10;  however,  this  damage  was  not  con- 
sidered serious. 

The  damage  to  the  test  items  on  Shot  10  was  considerably  more  severe 
than  had  been  anticipated..!/  The  upper  part  of  the  truss  section  ( 3.4a) 
failed  and  overturned.  The  overturning  began  when  the  vertical  posts 
and  diagonals  on  the  ground  zero  side  pulled  loose  from  the  gusset  plates. 
The  entire  upper  section  then  rotated  about  the  rear  posts,  which  even- 
tually failed.  The  rear  diagonals  remained  intact  even  after  the  upper 
section  reached  the  ground.  This  condition  is  shown  in  Figs.  4.1  through 

M- 

A small  permanent  set  was  observed  on  the  top  chord  component  3«^» 

The  majority  of  the  guy  cables  again  failed.  The  bottom  chord  of  struc- 
ture including  ties  and  rails,  was  undamaged.  The  bottom  chord 

(3.4c)  appeared  to  have  sustained  damage  from  flying  debris.  Damage 
was  confined  to  bending  of  the  gusset  plates  (see  Fig.  4.4). 

1/  The  test  was  not  originally  intended  for  Shot  10  precursor  conditions. 
When  it  was  known  that  the  test  structures  were  to  be  located  in  the 
precursor  region  of  Shot  10,  rough  response  computations  indicated  that 
some  damage  to  the  truss  section  might  occur,  but  not  to  the  degree 
actually  sustained. 
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Table  4.1  - Summary  of  Field  Conditions 


Shot 

No. 

Orientation 
(Degrees 
from  Normal) 

Over- 

Pressure 

(psi) 

* 

Duration 

(sec) 

Intended 

Mbs 

3.4a 

9 

2000 

2330 

20 

O.78 

1930 

4.6 

9.0 

0.57 

3.4b 

9 

2000 

2300 

21 

11.5 

0.78 

1925 

4.5 

9.0 

0.57 

3.4c 

9 

2000 

2275 

21 

11.6 

0.78 

1920 

4.5 

mm 

0.57 

3.4e 

9 

2000 

2230 

21 

0.77 

1915 

4.4 

9.0 

0.57 

3.4f 

9 

2000 

2200 

22 

n.8 

0.77 

10 

1910 

4.3 

9.0 

0.57 

Ambient  Preshot  conditions:  Shot  9 air  pre8SurePQ  = 13.2  psi; 

air  temperature  = l6.T°  C 
Shot  10  air  pressure  PQ  = 13-2  psi; 
air  temperature  = 14.8°  C 

* Duration  of  the  positive  phase  of  the  blast  wave 

Data  obtained  from.  Summary  Report  of  the  Technical  Director. 
WT-782. 
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4.3  STRAIN  GAGE  DATA 


All  of  the  strain  gages  provided  records  of  some  sort  on  both 
Shots  9 and  10.  However,  five  records  on  Shot  9 and  all  eighteen  re- 
cords on  Shot  10  are  considered  to  be  unusable.  Comments  on  the  strain 
gage  data  for  Shot  9 are  shown  in  Table  4.2  and  discussed  below.  Typi- 
cal records  are  shown  in  Fig.  4.5.  All  records  on  Shot  10  consisted 
entirely  of  what  might  be  termed  high-frequency  hash;  many  of  these 
show  predominant  zero  shifts  and  off-scale  readings.  It  is  believed 
that  electromagnetic  effects  associated  with  the  detonation  of  the  bomb 
were  responsible  for  the  erratic  gage  behavior  in  Shot  10.  This  be- 
havior is  considered  in  the  UPSHOT-KNOTHOLE  Project  3.28.1,  Final  Re- 
port ("Structures  Instrumentation,"  WT-738) . 

The  records  were  played  back  by  BRL  at  two  speeds;  a slow  speed 
which  gave  information  during  and,  in  most  cases,  well  beyond  the  posi- 
tive loading  period,  and  a fast  speed  which  indicated  greater  detail  of 
the  response  during  the  first  300  ms  or  so  after  shock  arrival.  Linear! 
zatlons  were  made  of  both  playbacks.  The  entries  in  Table  4.2  labeled 
ECq  and  EC2  refer  to  electrical  calibration  marks;  EC]_  is  a field  cali- 
bration taken  just  prior  to  shot  time,  and  SC2  is  an  earlier  calibration 
taken  some  time  before  the  tests.  No  postshot  calibrations  were  made. 
As  is  seen,  these  marks  are  generally  only  approximately  equal,  and  in 
those  cases  marked  "not  equal,"  the  difference  is  of  the  order  of  the 
strain  calibration  jumps.  The  linearizations  were  based  on  BRL's  best 
estimate  of  the  proper  calibration;  where  ECp  differed  appreciably  from 
ECg,  the  possible  error  involved  was  not  estimated  by  BRL. 

Inspection  of  the  Shot  9 records  points  out  the  strong  possibility 
that  at  least  three  significant  baseline  shifts  occurred  in  all  records 
between  about  40  and  100  ms  after  shock  arrival.  All  of  the  strain 
channels  utilized  the  same  recorder  unit.  Tracings  of  the  original 
fast  playbacks  from  structures  3* 4a  and  b are  shown  in  Fig.  4.6.  What 
appear  to  be  three  baseline  shifts  are  seen  to  occur  at  the  same  ab- 
solute time  (i.e.,  time  measured  from  the  bomb  detonation).  The  ap- 
proximate time  increment  between  these  shifts  for  the  Shot  9 records 
are  listed  in  Table  4.2.  The  average  time  between  jumps  is  about  10 
and  55  ms;  the  uncertainty  in  establishing  these  times  from  inspection 
of  the  records  is  of  the  order  of  variation  in  the  values  listed  in 
Table  4.2. 

On  the  average  these  jumps  occur  at  1786,  1796,  and  1851  ms  after 
bomb  detonation.  As  indicated  in  the  table,  of  the  apparently  usable 
records  only  record  3.4eS2  does  not  show  three  jumps,  but  this  record 
is  suspect  because  of  the  extremely  low  signal  level.  On  some  records 
(e.g.,  3.4aS-[_)  it  is  not  at  first  obvious  that  the  jumps  occurring  at 
the  above  times  are  indeed  baseline  shifts;  only  when  all  records  are 
viewed  in  this  light  does  this  conclusion  seem  inescapable,  see  Fig.  4.6 

The  magnitude  of  the  first  baseline  jump  is,  on  the  average,  about 
+80  per  cent  of  the  maximum  recorded  strain;  for  the  second  jump  it  is 
about  -40  per  cent  of  the  maximum  strain;  and  for  the  third  jump  it  is 
also  about  -40  per  cent.  The  fact  that,  on  the  average,  these  jumps 
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Table  4.2  - Summary  of  Strain  Gage  Data,  Shot  9 


Strain 

Record 

Remarks 

rime  Increment 
Between  Probable 
Baseline  Shifts 
(ms) 

3.4aS! 

ecl  i ec2 

T 

3.4aS2 

EC]_  i EC2 

3.4aS3 

EC^  - EC2 

3.4aS^ 

EC]_  = EC2 

ESQ 

3.4bSL 

EC1  - EC2 

11,  55 

3-4bS2 

ECX  i EC2 

10,  55 

3.4bS 

Extremely  hashy  baseline;  signal  off 

3 

scale,  not  usable 

3.4bS4 

EC1^  EC2 

12,  54 

3.4cS1 

EC1  — ECg 

11,  54 

3.4cS2 

Zero  shift;  no  recognizable  signal 

3.4cS3 

EC^  EC2 

11,  54 

3.4cS^ 

EC1'^  EC2 

10,  54 

3.4eSx 

High  frequency  hash,  not  usable 

3.4eS2 

EC1  = EC2,  Extremely  low  signal  level 

3.4fSx 

High  frequency  hash,  not  usable 

3*4fS2 

Zero  shift;  no  recognizable  signal 

3.4fS3 

EC1  - ECg 

3.4fS4 

ECX  ^ ECg 

99H 

Tig.  4.4  Postshot,  Bending  of  Gusset  Plate,  Bottom  Chord,  3*4e 


Strain  €,  in. /in.  Strain  £,  fJ.  in. /in.  Strain 


0 0.8  1.6  2.4  3.2  4.0  4.8  5*6  6.4 

Time,  t,  (sec) 

a.  3.4aSl,  Shot  9 (Slew  Playback) 


0 40  80  120  160  200  240  280  320  360  400  440  480  520  560 

Tine,  t,  (ns) 

b.  3.4aSl,  Shot  9 (Past  Playback) 


40  80  120  160  200  240  280  320  360  400  440  480 

Tine,  t,  (ns) 

c.  3.4b8l,  Shot  9 (Fast  Playback) 

Fig.  4.5  Typical  Strain  Gage  Records 


seem  to  nullify  each  other  suggests  the  possibility  that  the  records 
can  be  corrected  in  a simple  fashion  prior  to  the  time  of  the  third 
jump,  and  may  be  acceptable  or  they  stand  beyond  this  time.  This  was 
done  in  the  analysis  of  certain  of  the  records,  as  discussed  in  Chap- 
ter 5.  However,  the  origin  of  these  jumps  is  not  understood,  and  their 
apparent  nullification  may  be  entirely  coincidental. 

The  strain  calibration  steps  are  indicated  on  the  original  BRL 
playbacks  (see  Fig.  4.6).  A spot  check  of  the  linearized  data  against 
this  calibration  has  indicated  certain  inconsistencies  which  ARF  has 
not  been  able  to  resolve,  and  which  seriously  affect  the  analysis  of  at 
least  one  strain  record. 

In  particular,  the  linearized  record  for  strain  gage  3.4aSi,  Shot  9, 
does  not  agree  with  the  calibration  steps  given  on  the  original  play- 
back of  this  gage;  the  maximum  discrepancy  is  of  the  order  of  40  per 
cent.  Figure  4.7  compares  the  calibration  curve  for  this  record  as 
determined  from  the  calibration  steps  on  the  original  playback,  and  the 
calibration  curve  which  was  undoubtedly  used  in  the  BRL  linearization 
procedure.  This  latter  curve  was  determined  by  comparing  obvious  time 
details  between  the  BRL  linearization  and  the  original  playback.  It  is 
indeed  curious  that  the  apparent  calibration  curve  is,  as  well  as  can  be 
told,  a straight  line  passing  through  zero  and  the  maximum  calibration 
step,  whereas  the  curve  obtained  directly  from  the  calibration  steps  shows 
a characteristic  non-linear  shape. 

The  3*4aS1  record  was  replotted  according  to  this  latter  curve  and 
the  subsequent  analysis  was  carried  out  for  both  linearizations.  As 
discussed  in  Section  5,1.2,  the  fact  that  the  revised  ARF  linearization 
compares  favorably  with  the  other  records  from  the  3.4a  structure  lends 
some  support  to  the  correctness  of  this  revision.  Discrepancies  of  this 
type  were  noted  in  other  gage  records  but  the  percentage  deviation  did 
not  exceed  about  15  per  cent. 


48 


iMcussm 


CHAPTER  5 


POST-TEST  CONSIDERATIONS 


5.1  INTERPRETATION  OF  THE  STRAIN  DATA 


5.1.1  Introduction 


Visual  inspection  of  the  strain  data  did  not  yield  any  obvious 
information  as  to  the  nature  of  the  applied  blast  forces.  As  discussed 
in  Section  3*3.2,  this  vas  not  unanticipated  and  the  numerical  data  re- 
duction scheme  referred  to  as  transient  analysis  vas  applied  to  all 
usable  strain  records.  This  work  is  discussed  in  Section  5.1.2. 

Due  to  the  relatively  unsuccessful  application  of  transient 
analysis,  another  approach  to  the  data  reduction  problem  vas  employed 
in  an  attempt  to  get  additional  useful  information  from  the  data.  This 
method,  referred  to  as  graphical  analysis,  is  presented  in  Section  5.1.3. 

UPSHOT-KNOTHOLE  Project  3*3  ("Tests  on  the  Loading  of  Horizontal 
Cylindrical  Shapes,"  WT-722)  vas  concerned  with  exactly  the  same  data 
reduction  problems  as  in  the  present  program.  Toward  the  end  of  the 
post-test  phase  of  Project  3*3,  work  was  initiated  on  a rapid  trial- 
and-error  procedure  of  curve  fitting  utilizing  an  analog  computer. 

While  this  technique  was  not  applied  to  the  present  Project  3.4  strain 
data,  a brief  description  of  the  method  is  contained  in  Section  5.1.4. 

5.1.2  Application  of  Transient  Analysis 

It  was  originally  felt  that,  due  to  the  relatively  complicated 
response  of  the  dynamic  systems  involved,  the  magnitude  and  time  details 
of  the  incident  blast  loading  probably  could  not  be  determined  adequately 
from  inspection  of  the  strain  records  alone.  Furthermore,  it  did  not 
appear  feasible  to  determine  the  response  of  most  of  the  test  items  to 
even  a known  input  by  analytical  procedures.  Thus,  the  analysis  of  the 
data  could  not  proceed  by  ordinary  means.  This  situation  was  realized 
during  the  pretest  planning  phase  of  the  program  and  a numerical  data- 
reductlon  scheme,  referred  to  as  "transient  analysis,"  was  developed 
which  held  the  promise  of  determining  the  forcing  function  from  a known 
response  even  when  the  transfer  functions!/  of  the  system  were  not  known. 
1 / The  transfer  function  refers  to  the  response  of  the  system  in  one 
mode  due  to  a unit  Impulse  applied  in  another  mode  ( see  Appendix  A) . 
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In  general,  the  solution  to  this  problem  is  not  unique,  even  for 

a-lifear  dynamic  system>  tbat  is>  a given  response  in  a particular  mode 
of  the  system  may  result  from  more  than  one  combination  of  forcing  func- 
tions in  the  various  modes.  At  the  outset  it  was  believed  that,  in  the 
present  application,  uniqueness  stemmed  from  the  fact  that  both  the  un- 
known  forcing  and  transfer  functions  could  be  adequately  represented  by 
a particular  class  of  analytic  functions.  Specifically,  the  forcing 
„u^c_^on  was  represented  by  a finite  sum  of  exponential  terms  (i.e., 
e functions; , all  of  whose  exponents  and  coefficients  were  real:  the 
transfer  function  was  represented  by  a similar  series  all  of  whose  ex- 
ponents and  coefficients  occurred  in  complex  conjugate  pairs.  (This  is 
equivalent  to  assuming  a forcing  function  which  has  no  harmonic  compo- 
nents, and  transfer  functions  consisting  of  damped  sinusoids.)  Based 
on  this  particular  representation,  the  method,  in  brief,  consisted  of 
curve-fitting  the  known  response  (i.e.,  strain  data)  to  a similar  ex- 
ponential series  and  then  identifying  those  coefficients  and  exponents 
belonging  to  the  forcing  function.  In  the  process  the  transfer  func- 
tion is  also  determined.  The  detailed  development  of  the  method  of 
transient  analysis  as  well  as  a critical  discussion  of  its  general  appli- 
cability is  presented  in  Appendix  A.  The  present  section  is  concerned 
with  the  specific  application  of  the  method  to  the  strain  data. 

The  analysis  was  applied  to  all  four  records  on  structure  3.4a 
records  S^,  §2  s4  °n  and  records  Si,  So  and  Sk  on  3.4c  - all 

on  Shot  9.2/.  These  strain  records  were  fit  to  a five- term  exponential 
series  as  described  in  Appendix  A.  It  was  hoped  that  the  fifth-order 
polynomial  resulting  from  this  fit  would  have  three  real  roots  and  two 
complex  conjugate  roots.  If  so,  the  analysis  would  yield  a forcing 
function  consisting  of  three  exponential  terms  and  a transfer  function 
containing  a single  frequency.  However,  in  all  but  one  case  the  analysis 
provided  only  a one  term  approximation  to  the  forcing  function  and  two 
frequencies.  The  results  obtained  are  shown  in  Table  5.1  with  only  a 
relative  scale  being  indicated  for  the  magnitude  of  the  forcing  function 
since  a true  force  scale  cannot  be  determined  from  the  analysis  alone 
(See  Appendix  A.) 

The  one  case  in  which  a three- term  approximation  to  the  forcing 
function  was  obtained  (i.e.,  for  record  3.4aS2)  is  plotted  in  Fig.  5.1. 

The  shape  of  the  drag  curve  from  Fig.  5.12  (based  on  the  pretest  load 
prediction  scheme)  is  shown  for  comparison.  With  reference  to  Table  5.1, 
the  decay  rate  of  the  forcing  functions  determined  from  the  3.4a  strain 
records  appear  to  be  reasonable  in  terms  of  what  is  known  about  the  load- 
ing. Inspection  of  these  strain  records  shows  a frequency  of  about  29  cps 
and  possibly  one  of  about  50  cps,  (see  Fig.  4.5)  which  is  of  the  order  of 
magnitude  found  for  the  frequencies  of  the  transfer  function.  That  the 
analysis  completely  missed  the  principal  frequency  of  about  1.5  cps. 


2/  The  analysis  was  applied  to  these  records  before  the  probable  oc- 
curence of  baseline  shifts  and  calibration  errors  was  appreciated 
(see  Section  4.3).  The  possible  error  introduced  in  the  analysis 
is  considered  later  in  this  section. 
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however,  is  not  surprising  since  the  results  are  based  on  data  for  only 
the  first  300  to  400  ms  of  the  loading.  The  results  for  the  other  struc- 
tures shown  in  Table  5.1  indicate  that  the  proper  frequencies  may  have 
been  picked  out,  but  most  of  the  decay  rates  appear  to  be  too  large. 

In  general,  the  results  indicated  in  Table  5.1  are  not  believed 
to  be  an  adequate  interpretation  of  the  strain  data  and,  as  such,  it 
must  be  concluded  that  this  first  application  of  transient  analysis  was 
unsatisfactory.  In  an  attempt  to  understand  just  why  the  analysis  failed 
a rather  detailed  post-test  reappraisal  of  the  method  was  undertaken. 

Much  of  the  critical  discussion  presented  in  Appendix  A is  a result  of 
this  additional  study. 

The  first  application  of  transient  analysis  was  completed  before 
the  occurrence  of  probable  baseline  shifts  and  calibration  errors  in 
the  strain  data  was  fully  appreciated  (see  Section  4-. 3) . The  fact  that 
the  strain  scale  of  the  records  might  be  in  error  by  a constant  factor 
would  not  affect  the  reliability  of  the  analysis,  inasmuch  as  the  forc- 
ing function  can  be  determined  only  within  an  arbitrary  scale  factor. 
However,  if  the  time  details  of  the  strain  record  are  distorted  due  to 
non-linear  scale  errors  or  baseline  shifts,  the  situation  is  different. 
Here  the  forcing  function  determined  by  transient  analysis  might  not  have 
physical  meaning  since  it  would  be  based  on  prime  data  that  was  in  part 
arbitrary. 

Thus,  the  uncertainty  attributed  to  the  strain  data  could  account 
to  some  extent  for  the  unsatisfactory  results  obtained,  although  present 
understanding  of  the  transient  analysis  makes  it  impossible  to  determine 
the  extent  to  which  errors  in  the  prime  data  can  influence  the  final 
results.  In  the  present  instance,  however,  it  is  felt  that  the  failure 
of  this  first  application  of  transient  analysis  can  be  more  readily  ex- 
plained in  terms  of  the  poor  curve-fit  obtained  for  the  known  response. 
Whether  it  was  due  to  computational  errors  or  simply  an  inadequate  num- 
ber of  terms  in  the  approximating  series  (most  probably  the  latter) , 
the  fact  remains  that  the  determination  of  the  forcing  functions  was 
based  on  analytical  expressions  which  were  not  adequate  representations 
of  the  measured  strain  time  curves  - even  though  these  curves  may  them- 
selves be  in  error. 

As  discussed  in  Appendix  A,  the  valid  application  of  transient 
analysis  is  limited  to  linear  systems  where  the  forcing  functions  as- 
sociated with  each  mode  are  either  zero  or  timewise  proportional.  When 
the  forcing  functions  are  not  all  timewise  proportional,  the  results 
of  the  analysis  probably  cannot  be  given  physical  meaning.  (This  limi- 
tation of  the  method  was  not  fully  appreciated  at  the  outset.)  Now  it 
is  possible  that  in  as  complicated  a structure  as  the  truss  bridge,  the 
forcing  functions  associated  with  each  mode  are  not  timewise  proportional 
But  this  cannot  be  demonstrated  one  way  or  the  other  at  the  present  time, 
and  the  applicability  of  the  method  is  not  sufficiently  well-understood 
to  decide  definitely  as  to  its  validity. 

A second  application  of  the  method  of  transient  analysis  based 
on  a more  realistic  fit  of  the  strain  records  was  initiated,  but  this 
work  could  not  be  completed  within  the  scope  of  the  post-test  phase  of 
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the  program.  Only  preliminary  results  with  respect  to  one  strain  re- 
cord were  obtained,  and  while  these  appear  encouraging  no  general  con- 
clusions can  be  drawn.  The  data  considered  was  strain  record  3.,4&S^  on 
Shot  9«  A comparison  is  made  between  the  actual  strain  record!/  and 
the  results  of  the  curve-fitting  in  Fig.  5*2.  This  rerun  of  the  analy- 
sis differed  from  the  first  computation  in  that  the  approximating  series 
consisted  of  seven  exponential  terms  rather  than  five,  and  the  obvious 
frequency  of  about  1.5  cps  was  cui ve-f itted  directly.  The  second  fre- 
quency as  determined  by  the  analysis  itself  was  2 9 cps,  which  agrees 
very  closely  with  an  observable  frequency  in  the  data. 

While  the  results  of  this  curve-fitting  appear  to  be  realistic 
(see  Fig.  5.2),  they  could  no  doubt  be  improved.  The  resulting  forcing 
function  is  shown  in  Fig.  5*7  and  is  compared  with  certain  results  ob- 
tained in  the  following  section.  At  present,  the  curve  labeled  "graphi- 
cal analysis"  in  Fig.  5.7  may  be  accepted  as  representing  an  independent 
estimate  of  the  forcing  function  component  leading  to  the  measured  strain 
of  gage  3.4aSi  (see  Section  5.1.2).  Inasmuch  as  the  forcing  function 
obtained  from  transient  analysis  has  an  arbitrary  scale  factor,  the  pre- 
sent results  were  fitted  to  the  results  of  graphical  analysis  at  about 
100  ms. 

With  reference  to  Fig.  5-7,  the  forcing  function  determined  by 
transient  analysis  shows  the  same  broad  time  details  as  does  the  graphi- 
cal approach.  That  is,  the  time  of  rise  to  maximum  load  and  the  rate 
of  decay  during  the  drag  phase  of  the  loading  are  in  good  agreement. 
However,  the  time  details  during  the  diffraction  period  are  not  adequately 
reproduced  and,  indeed,  such  detail  could  not  be  expected  from  a three- 
term  exponential  series. 

On  the  basis  of  this  comparison  transient  analysis  appears  to  be 
a feasible  method  in  the  present  application,  provided  that  the  analy- 
tical approximation  to  the  strain  data  (i.e.,  curve-fitting)  is  itself 
realistic.  But  this  is  not  a definite  conclusion,  and  much  additional 
work  remains  to  be  done  in  application  of  the  method. 

5.I.3  Graphical  Analysis 

In  view  of  the  rather  unsuccessful  results  of  the  first  application 
of  transient  analysis,  it  was  necessary  to  follow  smother  approach  in  an 
effort  to  obtain  the  desired  information  from  the  strain  data.  The  meth- 
od adopted  was  an  attempt  to  remove  the  obvious  harmonic  components  of 
the  response  by  inspection  and  thereby  yield  at  least  the  shape  of  the 
forcing  function.  For  lack  of  a better  name  this  technique  was  referred 
to  as  graphical  analysis.  The  method  was  applied  with  greatest  success 
to  the  3.4a  strain  records,  and  only  these  results  are  considered  in  this 
section. 

The  probable  baseline  shifts  in  the  3.4a  strain  records  (see  Sec- 
tion 4.3)  were  eliminated  prior  to  the  application  of  the  present  analysis. 
y The  baseline  shifts  on  this  particular  record  were  not  appreciable 
and,  hence,  they  were  not  compensated  for  (see  Fig.  4.6). 
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The  procedure  consisted  of  determining  both  the  time  of  occurrence  and 
magnitude  of  the  shifts  by  visual  inspection  of  the  linearized  strain 
records,  and  replotting  the  records  accordingly.  Three  baseline  shifts 
were  removed  from  each  of  the  four  strain  records  in  this  fashion. 

The  strain  records  from  the  3.4a  truss  bridge  section  show  a well- 
defined  frequency  of  about  1.5  to  1.6  cps  (Fig.  4.5)  which  is  undoubted- 
ly due  to  vibration  of  the  top  chord  members  relative  to  the  bottom  chord 
(see  Section  5*2).  The  period  of  this  vibration  is  Just  slightly  less 
than  the  duration  of  the  positive  phase  of  the  blast  wave. 

It  was  possible  to  determine  an  analytical  approximation  to  this 
principal  harmonic  component  of  the  response  by  means  of  measurements 
taken  from  the  strain  records.  The  analytical  approximation  determined 
from  each  of  the  strain  records  is; 

for  strain  record  3.4aSi 

l49e"°*3lt  cos  (9-96t  - 1.20) 


97e'°*31t  cos  (9.96t  - 1.20) 

for  strain  record  3.4aS2  ^ (5*1) 

97e"°*31tcos  (9.96t  - 1.20) 
and  for  strain  records  3«4aS2  and  3.4aS4 

97e"°*31t  cos  (9«96t  + 2.29),  t in  sec. 

The  two  amplitudes  given  for  strain  record  3.4aS^  are  based  on  the  BRL 
linearization  (amp  = 149)  and  the  revised  ARF  linearization  (amp  = 97), 
as  discussed  in  Section  4.3.  The  numerical  values  given  above  are  all 
average  values,  the  variation  in  each  being  about  + 15  per  cent.  The 
fact  that  the  four  records  yield  the  same  amplitude  if  the  revised 
linearization  is  accepted,  lends  some  support  to  this  procedure. 

The  first  step  in  the  analysis  was  to  subtract  (at  each  instant 
of  time)  the  harmonic  component  given  by  Equation  5.1  from  the  appro- 
priate strain  record.  The  results  of  this  subtraction  are  referred  to  in 
Figs.  5.3  through  5.6  as  the  "first  reduction."  From  these  results  a 
rather  obvious  median  curve  can  be  faired  in,  about  which  the  record 
appears  to  be  oscillating.  This  curve  is  referred  to  as  the  "second 
reduction"  in  Figs.  5*3  through  5*6  and  represents,  according  to  this 
approach,  the  shape  of  the  forcing  function  contributing  to  the  reaction 
in  the  individual  sensor  bar. 

In  the  case  of  record  3.4aS]_  a second  frequency  can  be  identified 
in  both  the  BRL  and  ARF  linearizations,  which  permits  a somewnat  more 
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!°,nthe  second  reduction.  From  Fig.  5.3  this  frequency 
is  seen  to  be  about  29  cps  and,  as  discussed  in  Section  5.2,  is  due  ao- 

i?  t?e  VibraJi0n  °f  the  entire  &tructure  the  sensor  supports. 
AflF  r'SWSe  tOT  «“ 

132e"7-6t  cos  l82t 

This  component  was  subtracted  from  the  first  reduction  and  the  resulting 
second  reduction  for  this  record  is  shown  in  Fig.  5.7  where  the  strain^ 
scale  has  been  converted  to  force  per  unit  frontal  area  of  the  upstream 
^rUS!C  reduction  introduces  a strong  harmonic  component  dur- 

ing  the  first  50  ms  of  the  loading  and  then  becomes  increasingly  similar 
to  the  first  reduction  as  the  29  cps  frequency  is  damped  out. 

,.The  Predicted  diffraction  loading  from  Fig.  5.12,  is  compared 
with  the  result  of  the  second  reduction  in  Fig.  5.7.  As  can  be  seen, 
the  two  curves  compare  quite  favorably  when  allowance  is  made  for  the 
uncertainty  in  zero  times.  Thus,  these  results  tend  to  substantiate 
the  load  prediction  scheme  during  the  diffraction  period,  at  least  with 
respect  to  time  details. 

As  discussed  in  the  previous  section,  the  rerun  of  transient  analy- 
sis was  carried  out  for  strain  record  3.4083.,  Shot  9.  These  results 
are  compared  with  the  results  of  the  present  method  in  Fig.  5.7.  Inas- 
much  as  the  forcing  function  obtained  from  transient  analysis  has  an 
arbitrary  scale  factor,  the  solid  curve  in  Fig.  5.7  was  fitted  to  the 
present  results  at  about  100  ms.  A discussion  of  this  comparison  is 
contained  in  the  previous  section. 

The  horizontal  component  of  force  transmitted  to  the  sensor  bars 
is  proportional  to  the  algebraic  difference  in  strain  as  given  by  this 
equation.  According  to  this  equation,  which  results  from  treating  the 
structure- sensor  system  as  a statically  determinate  problem,  the  hori- 
zontal  force  may  be  obtained  from  either  strain  records  Si  and  Sp  or 
records  S3  and  S4.  The  horizontal  force  per  unit  projected  area  of  the 
first  truss  as  determined  from  strain  records  S3.  and  S2  is  shown  in 
Fig.  5.8.  Results  are  shown  for  both  the  BRL  and  ARF  calibration  of 
strain  record  SL.  The  horizontal  force,  as  determined  from  records  So 
and  S4,  is  shown  in  Fig.  5.9.  3 

The  horizontal  force  as  determined  from  strain  records  Si  and  So 
is  compared  with  the  predicted  horizontal  loading  based  on  the  measured 
blast  parameters  for  Shot  9 (see  Section  5.2)  in  Fig.  5.8.  Allowing  for 
uncertainties  in  the  location  of  zero  time,  these  results  utilizing  the 
ARF  linearization  of  record  S3.  are  seen  to  compare  remarkably  well  with 
pretest  predictions;  the  results  using  the  BRL  calibration  indicate  the 
same  time  details,  but  show  pressures  about  twice  the  predicted  values 

— /'  doping  term  is  admittedly  quite  approximate,  since  only  several 

cycles  of  the  vibration  can  be  identified  in  Fig.  5.4. 
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during  the  drag  phase.  The  horizontal  force  computed  from  records  S3 
and  S4  (see  Fig.  5.9)  shows  different  time  details  during  the  diffrac- 
tion period  and  much  larger  drag  forces.  In  this  respect  the  drag 
forces  are  in  agreement  with  the  results  of  the  BRL  linearization  of 
record  Si  shown  in  Fig.  5*8.  The  differences  in  horizontal  force  as 
determined  from  the  front  and  rear  sensor  gages  may  be  due  in  part  to 
orientation  effects.  The  predicted  loading  is,  of  course,  based  on 
head-on  orientation;  wheras  the  blast  wave  which  struck  the  3.4a  bridge 
section  was  oriented  about  20  degrees  from  the  normal. 

In  view  of  the  uncertainty  attending  both  the  pretest  load  pre- 
diction scheme  and  the  strain  data,  the  agreement  between  measurement 
and  shock  theory  in  Fig.  5*8  is  indeed  remarkable,  and  may  in  fact  be 
coincidental.  Considering  the  rather  conflicting  nature  of  other  re- 
sults obtained  however,  one  is  not  justified  in  definitely  concluding 
that  the  test  data  support  the  pretest  load  prediction  method. 

5.1.4  Analog  Computer  Method 

UPSHOT-KNOTHOLE  Project  3.3  ("Tests  on  the  Loading  of  Horizontal 
Cylindrical  Shapes','  WT-722)  was  concerned  with  many  of  the  same  problems 
encountered  here.  That  is,  the  test  was  designed  to  determine  the  load- 
ing on  horizontal  cylinders  utilizing,  in  part,  a net  force  measurement 
system  with  strain  output  identical  to  the  one  in  the  present  test.  In 
connection  with  Project  3.3,  work  was  initiated  on  the  use  of  the  ARF 
analog  computer  to  deduce  the  desired  information  from  the  strain  re- 
cords. While  time  would  not  allow  this  technique  to  be  employed  in  the 
present  program,  this  approach  appears  to  be  well  suited  to  it  and  will 
be  summarized  briefly  in  this  section. 

The  method  of  solution  utilizing  the  analog  computer  consists  of 
simulating  a dynamic  system  on  the  computer  which  yields  displacements 
similar  to  the  strain  record  under  consideration,  and  then  adjusting  the 
circuit  parameters  so  that  the  two  displacement-time  curves  are  essen- 
tially the  same.  The  system  equations  are  taken  to  be  a series  of  linear 
second-order  differential  equations  with  constant  coefficients,  an  as- 
sumption which  is  also  basic  to  the  method  of  transient  analysis.  The 
forcing  function  may  be  represented  by  a variety  of  analytical  or  graphi- 
cal functions  on  the  computer  but  only  decaying  exponential  terms  were 
utilized  in  connection  with  the  work  on  Project  3*3>  again  as  in  the 
method  of  transient  analysis. 

On  first  thought,  this  approach  might  seem  rather  hopeless  in 
view  of  the  large  number  of  unknown  quantities  involved,  i.e.,  the  con- 
stants of  the  system  equations  and  the  parameters  of  the  forcing  func- 
tion. However,  in  practice,  the  method  appears  to  be  quite  feasible, 
and  while  not  carried  sufficiently  far  as  to  provide  any  conclusive 
results,  gives  indications  of  being  an  advantageous  approach.  The  fea- 
sibility of  the  method  in  practice  stems  from  the  fact  that  the  strain 
records  to  be  matched  can  be  plotted  on  graph  paper  and  the  displacement 
arising  from  the  simulated  system  can  be  superposed  automatically  be 
means  of  an  "x-y  recorder."  This  allows  a visual  estimate  of  the  ac- 
curacy of  the  matching,  and  permits  much  more  judgment  on  the  part  of 
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the  operator  in  carrying  out  the  solution  than  is  possible  when  using 
a numerical  method.  This  advantage,  together  with  the  brief  time  re- 
quired for  varying  parameters  and  obtaining  solutions,  permits  the 
matching  to  be  completed  in  a relatively  short  time. 

Use  of  the  analog  computer  technique  as  a means  of  data  reduc- 
tion is  only  in  the  very  early  preliminary  stages.  However,  this  work 
has  indicated  a reasonable  expectation  of  ultimate  success,  and  it  is 
strongly  recommended  that  this  work  be  pursued  further  as  a part  of  an 
over-all  study  into  the  data-reduction  problems  associated  with  net  force 
measurement  systems.  One  advantage  of  the  analog  computer  approach  which 
concerns  the  forcing  function  approximation  should  be  emphasized.  It 
is  not  necessary  in  the  analog  scheme  to  be  restricted  to  an  exponential 
approximation  of  the  forcing  function  since  it  is  possible  to  generate 
many  other  functions  on  the  computer.  For  example,  one  form  of  repre- 
sentation which  could  be  handled  easily  is  a series  of  straightline 
sections  during  the  initial  loading  period,  followed  by  either  a linear 
or  exponential  decay  to  zero.  Existing  load  prediction  schemes  could 
be  generated  on  the  computer  and  the  response  of  the  simulated  system 
compared  with  the  actual  response  data. 

5.2  RESPONSE  OF  THE  TRUSS  BRIDGE  SECTION,  3.4a 


5.2.1  Introduction 

The  truss  bridge  section  (3.4a)  suffered  a small  permanent  set 
of  about  3 in.  on  Shot  9 (no  precise  measurement  of  this  deformation 
was  taken),  and  was  completely  overturned  as  a result  of  Shot  10.  An 
attempt  is  made  in  this  section  to  correlate  the  predicted  blast  load- 
ing with  the  observed  response  of  the  bridge. 

The  truss  bridge  represents  a rather  complicated  structure,  and 
a precise  determination  of  its  response  to  even  a known  loading  is  not 
practical  on  the  basis  of  present-day  knowledge.  However,  circumstances 
are  such  that  an  approximate,  though  fairly  realistic,  estimate  of  the 
response  can  be  made  as  a result  of  the  test  data  obtained. 

The  strain  records  on  Shot  9 indicate  the  bridge  to  be  responding 
in  two  principal  modes;  one  with  a frequency  of  about  1.5  cps  which  is 
undoubtedly  due  to  vibration  of  the  top  chord  members  relative  to  the 
bottom  chord,  and  one  with  a frequency  of  abo'.t  29  cps  which  appears  to 
be  due  to  the  total  mass  of  the  bridge  vibrating  on  the  flexible  sensor 
supports.  In  the  following  discussion  the  higher  frequency  mode  will 
be  ignored  and  the  bridge  assumed  to  be  a single-degree-of-freedom  system 
possessing  a frequency  of  1.54  cps  of  a period  of  0.65  spc.£/  The  meas- 
ured frequency  determines  the  ratio  of  effective  stiffness  to  effective 
mass  of  the  bridge  section.  These  quantities  are  computed  in  the  fol- 
lowing section. 

£/  This  is  an  average  value,  the  variation  being  about  + 0.02  sec 
cycle-to-cycle  and  record- to-record. 
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5.2.2  Effective  Mass  and  Resistance  of  Bridge 


Figure  4.1  is  a postshot  10  photograph  of  the  bridge  showing  the 
manner  of  failure.  It  is  seen  that  the  members  above  the  lower  sway 
brace  strut  remained  essentially  rigid  and  that  plastic  hinges  formed 
at  the  ends  of  the  vertical  posts  just  below  this  strut.  It  seems  rea- 
sonable, therefore,  to  take  the  effective  mass  of  the  bridge  as  being 
the  total  mass  of  the  members  above  the  lower  sway  brace  strut  plus 
1/4  the  mass£/  of  the  unsupported  posts  and  main  diagonals.  This  mass 
is  found  to  be,  m = 142  lb-sec^/in.  Using  the  following  relationship 
for  the  period  of  a simple  mass-spring  system. 


T = 2 7T 


= 0.65  spc 


The  effective  stiffness  is  computed  to  be  ki  = 13.2  x 10^  kips/in. 

There  is  now  the  problem  of  determining  the  effective  resistance 
of  the  bridge  from  this  stiffness  value.  The  major  uncertainty  concerns 
the  effectiveness  of  the  main  diagonals  which  are  the  same  section  as 
the  posts  (14WF78)  but  of  longer  length,  and  oriented  to  bend  in  the 
strong  direction.  If  the  diagonals  are  discounted  for  the  moment,  and 
the  posts  assumed  to  act  as  cantilever  beams  with  guided  ends,  the  total 
stiffness  (for  four  posts)  is 


k - 48EI 

1 ‘ { 3 

„ 4 , „ K 

Using  I = 206  in.  for  a 14WF78  section  bending  in  the  weak  direction 
and  E = 30  x 10^  psi,  the  effective  length  X is  found  to  be  about 
23.5  ft.  This  is  close  to  the  unsupported  post  length  of  about  21  ft. 

In  computing  the  response  of  the  bridge,  the  resistance  function 
is  assumed  to  have  an  idealized  elastic-plastic  form  as  indicated  in 
Fig.  5.10.  The  stiffness  during  the  elastic  portion  is  given  by 


R 

k „ _R 
*1  X 

P 


where  Rp  is  the  maximum  (plastic)  resistance  and  Xp  is  the  limiting 
elastic  displacement.  Based  on  a yield  stress  of  O’y  = 45,000  psi,  the 
maximum  (plastic)  moment  of  a 14WF78  section  oriented  to  bend  in  the 
weak  direction  is  computed  to  be 


M 


<r  s 
y 


= 2.4 


x 10^, 


in. -lb  per  post. 


6/  This  is  the  usual  approximation  in  computing  the  fundamental  fre- 
quency of  a vibrating  cantilever  beam  with  a concentrated  mass  on 
the  end  wherein  the  effective  mass  of  the  system  is  taken  as  l/4 
(actually  0.23)  of  the  mass  of  the  beam  plus  the  total  end  mass. 
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For  a guided  cantilever  with  end  shear  R,  and  end  moment  M, 


R = (2 ji  ) M 

Using  this  relation  and  the  effective  post  length  computed  above,  the 
maximum  resistance  of  the  bridge,  R^,  is  computed  to  be 


„ 4 x 2 x 2.4  x 10°  . 

RP  • 68 

The  associated  yield  displacement  is  then 


X _ 68  x IQ3 

p kl  13-2  x 103 


= 5.2  in. 


The  maximum  plastic  resistance  determined  above  is  felt  to  be  a 
minimum  value,  inasmuch  as  the  resistance  of  the  diagonals  has  been 
neglected.  The  diagonals  are  also  14WF78  sections,  but  are  longer  (total 
length  about  50  ft)  and  are  oriented  to  bend  in  the  strong  direction. 

All  things  considered,  the  actual  resistance  is  not  believed  to  exceed 
the  computed  value  by  more  than  a factor  of  2,  and  this  range  of  resis- 
tance will  be  considered  in  the  following  response  computations. 

5.2.3  Response  Computations,  Shot  9 

Brooks  and  Newmark  ("Development  of  Procedures  for  Rapid  Compu- 
tation of  Dynamic  Structural  Response")  have  determined  the  response 
of  non-linear  mass-spring  systems  to  triangular- shaped  loading  functions 
which  are  typical  of  blast  loadings  based  on  current  prediction  methods. 
These  solutions  are  given  in  graphical  form  and  by  empirical  relation- 
ships determined  from  many  numerical  solutions.  In  particular,  Brooks 
and  Newmark  have  developed  an  empirical  formula  for  determining  the 
maximum  displacement  of  a mass-spring  system  having  an  elastic-plastic 
resistance  function  which  may  include  the  deadweight  overturning  effect 
of  the  structure,  see  Fig.  5.10.  The  loading  considered  consists  of  an 
initial  impulse  plus  an  initially  peaked  triangular  pulse.  This  rela- 
tionship is, 


A - D 

Xp  *2 


(5.2) 


where 


A = 2f-1  + * r-1 


and  where 


D = (27Tio/RpT)‘ 


F = peak  load  of  triangular  pulse 
m 

R = maximum  resistance 
P 

Xm  = maximum  displacement  of  mass 
Xp  = yield  displacement 

= slope  of  elastic  portion  of  resistance  function 
kg  = slope  of  plastic  portion  of  resistance  curve 

K1 

M = effective  mass  of  system 
T = 27r  = natural  period  of  system 


tg  = duration  of  triangular  load 
I = initial  impulsive  loading. 

Equation  5.2  becomes  exact  in  the  limiting  cases  for  both  in- 
finitely short  and  infinitely  long  duration  pulses,  i.e. , for  t2  > 0 
and  t2  > oo  . For  loads  of  an  intermediate  duration.  Brooks  and  New- 
mark  state  that  Equation  5.2  is  generally  accurate  to  within  6 per  cent. 
This  is  considered  completely  adequate  for  the  present  application. 

The  pretest  load  prediction  scheme  described  in  Chapter  3 was 
applied  to  the  3.4a  bridge  section  using  the  measured  blast  wave  para- 
meters for  Shot  9 (see  Table  4.1).  The  net  horizontal  load  is  shown 
in  Fig.  5.12.  Inspection  of  this  figure  shows  that  it  is  not  unreason- 
able to  approximate  the  loading  by  an  equivalent  initial  impulse  fol- 
lowed by  a triangular  decay.  The  shaded  area  under  the  net  load-time 
curve  from  zero  to  the  duration  of  the  diffraction  loading  (i.e.,  at 
t = 80  ms)  is  taken  to  be  the  initial  impulse  IQ  j/i  IQ  = 0.225  psi  sec._ 
The  remaining  portion  of  the  loading  (i.e.,  the  drag  loading)  is  ap- 
proximated by  a triangle  whose  initial  peak  (Fm)  is  taken  to  be  the 
peak  drag  force  and  whose  duration  (t]_)  is  chosen  so  that  the  total  drag 
impulse  is  preserved.  The  equivalent  loading  is  shown  in  Fig.  5*12; 

^m  = 3*2  Psi>  t2  = 0.32  sec. 

Equation  5.2  was  solved  numerically  for  the  maximum  resistance, 
Rp,  as  a function  of  maximum  displacement,  Xm,  based  on  a constant 


2/  Brooks  and  Newmark  state  that  a load  applied  over  a finite  time  in- 
terval may  be  approximated  by  an  equivalent  impulse  provided  the 
interval  is  less  than  about  l/lO  the  period  of  the  system.  In  the 
present  case  the  duration  of  the  diffraction  loading  is  about  12 
per  cent  of  the  period. 
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period,  T,  and  the  above  load  parameters.  The  results  are  plotted  in 
Fig.  5 • 13 . It  was  estimated  that  the  bridge  sustained  a 3-in.  permanent 
set  in  Shot  9,  or  a maximum  displacement  of  ^ = Xp  + 3 In*  From  Fig.  5*13> 
a displacement  of  this  amount  corresponds  to  a resistance  about  50  per 
cent  greater  than  the  estimated  minimum  value.  Since  the  maximum  pos- 
sible resistance  is  estimated  at  something  less  than  twice  the  minimum 
value,  this  result  is  plausible.  Based  on  the  estimated  minimum  value 
for  Rp,  Xm  = 15.4  in.  Xp  ,+  10  in.  For  the  estimated  maximum  value 
of  Rp,  Xm  = 9.7  in.  = XpX_/  - 0.7  in.  (i.e.,  the  bridge  would  have  re- 
mained elastic). 

A further  computation  was  carried  out  in  which  displacement  as 
a function  of  drag  coefficient  was  determined  for  the  range  of  interest 
in  resistance.  These  results  are  shown  in  Fig.  5.14,  where  the  drag 
coefficient  is  expressed  as  a percentage  of  the  nominal  value  used  in 
the  pretest  load  prediction  scheme.  This  nominal  value  is  = 3 = C§, 
based  on  the  horizontal  projected  area  of  the  first  truss  only  as  dis- 
cussed in  Chapter  3.  The  value  of  as  determined  from  the  condition 

that  Xm  = Xp  + 3 in.  is  seen  from  Fig.  5-14  to  vary  between  0.6  and  1.3 
for  the  range  of  resistance  considered. 

There  is  one  important  conclusion  to  be  drawn  from  the  above  anal- 
ysis; namely,  that  the  pretest  load  prediction  scheme  for  the  3*^a  truss 
bridge  in  conjunction  with  a simplified  response  analysis  leads  to  a 
valid  estimate  of  the  actual  behavior  of  this  structure  under  field  con- 
ditions. That  is  not  to  say  that  the  present  result  amount  to  a veri- 
fication of  the  load  prediction  scheme  or  lead,  say,  uo  any  revision  in 
over-all  drag  coefficient.  The  present  analysis  is  certainly  too  crude 
to  permit  this,  although  such  was  the  intent  of  uhis  test  program. 

Rather,  existing  load  and  response  prediction  methods  would,  in  the  pre- 
sent case,  lead  one  to  conclude  that  the  3*^a  bridge  would  not  suffer 
excessive  damage  for  Shot  9 conditions  — whether  the  displacement  turned 
out  to  be  5 in.  or  15  in.  — and  this  result  was  indeed  borne  out.  In 
other  respects  it  is  felt  that  this  analysis  has  no  quantitative  mean- 
ing. 


5.2.4  Response  Computations,  Shot  10 

The  3.4a  bridge  section  was  overturned  as  a result  of  the  loading 
experienced  in  Shot  10.  The  bridge  was  well  within  the  precursor  region 
in  this  shot  and  its  response  was  more  violent  than  had  been  anticipated, 
although  detailed  pretest  response  computations  had  not  been  made  for 
estimated  Shot  10  conditions.  One  result  of  the  UPSHOT-KNOTHOLE  Project 
3.1  ("Tests  on  the  Loading  of  Building  and  Equipment  Shapes,"  WT-721) 
was  to  propose,  quite  tentatively,  a load  prediction  scheme  applicable 
to  items  of  simple  geometry  in  the  precursor  region.  The  present  section 
is  concerned  with  response  computations  for  the  3.4a  bridge  based  on  this 
load  prediction  scheme. 

87  In  order  to  preserve  the  measured  period  of  O.65  spc,  the  yield 
displacement,  Xp  must  vary  directly  with  the  maximum  plastic  re- 
sistance, Rp.  For  Rp,],^  3 136  kip6,  Xp  = 10.4  in. 
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The  loading  method  developed  in  WT-721  stemmed  from  pressure 
gage  data  obtained  on  a solid  rectangular  parallelepiped  located  in 
the  precursor  region  of  Shot  10  that  rigidly  withstood  the  effects  of 
the  blast.  These  data  have  led  to  an  average  force-time  variation  ex- 
perienced by  the  structure  as  indicated  in  Fig.  5*15*  The  load  rises 
linearly  from  zero  to  a maximum  value  in  approximately  100  ms,  a time 
which  is  felt  to  be  independent  of  the  geometry  of  the  structure.  The 
load  then  decays  linearly  to  zero  in  a time  equal  to  the  positive  phase 
duration  of  the  blast  wave.  The  maximum  pressure  attained  is  taken  to 
be  the  peak  drag  pressure  associated  with  the  peak  overpressure  which 
would  have  occurred  over  an  "ideal"  surface  in  the  absence  of  a pre- 
cursor, multipled  by  some  over-all  drag  coefficient  for  the  structure. 

This  loading  is  markedly  different  from  what  is  generally  felt  to  the 
case  in  the  conventional  Mach  reflection  region  in  the  absence  of  a 
precursor.  There  is  a finite  rise  time,  no  diffraction  period,  and  a 
wave  form  during  the  drag  phase  which  differs  from  the  usually  accepted 
wave  form  of  the  dynamic  pressure  curve. 

There  exists  considerable  doubt  concerning  the  prediction  of  peak 
dynamic  pressure  in  the  precursor  region.  In  UPSHOT-KNOTHOLE  Project 
l.ld  ("Dynamic  Pressure  Versus  Time  and  Supporting  Air  Blast  Measurements," 
WT-71U),  dynamic  pressure  gages  were  located  on  the  main  blast  line  at 
a ground  range  of  1920  ft  and  elevations  of  10,  25,  and  40  ft  in  Shot  10. 
(The  3.4a  structure  was  at  a ground  range  of  about  1930  ft  and  extended 
from  15  ft  to  about  65  ft  above  ground  level.)  These  gages  recorded  max- 
imum dynamic  pressures  of  10.9,  11.6  and  12.4  psi,  in  order  of  increasing 
elevation.  However,  the  quantity  measured  by  these  gages  is  believed  to 
be  strongly  influenced  by  the  presence  of  suspended  dust  in  the  air  and, 
in  addition,  there  is  some  uncertainty  associated  with  the  calibration 
of  these  gages.  Theoretical  predictions  for  dynamic  pressure  over  an 
ideal  surface  are  contained  in  both  WT-714  and  WT-721,  and  indicate  a 
peak  pressure  of  about  9 psi  for  the  3.4a  bridge.  Using  a nominal  over- 
all drag  coefficient  of  3 for  the  bridge,  the  peak  force  per  unit  frontal 
area  of  the  first  truss  is  believed  to  be  in  excess  of  27  psi  as  indi- 
cated in  Fig.  5* 15* 

Brooks  and  Newmark  have  obtained  the  response  of  an  elastic-plastic 
mass-spring  system  to  the  triangular- shaped  loading  discussed  above  and 
shown  in  Fig.  5.15.  Their  results  are  presented  in  graphical  form  in 
"Development  of  Procedures  for  Rapid  Computation  of  Dynamic  Structural 
Response."  Based  on  this  solution,  the  response  of  the  3 -4a  truss 
bridge  section  has  been  computed  as  a function  of  the  peak  load  to  re- 
sistance ratio  and  is  shown  in  Fig.  5«l6.  Brooks  and  Newmark  have  re- 
sults for  k =-0.04  and  k = 0 (the  factor  k determines  the  relative 
shape  of  the  resistance  function,  see  Fig.  5.10  and  Equation  5.2).  In 
the  present  case,  k = -0.014  and,  as  seen  in  Fig.  5.1 6,  there  is  some 
uncertainty  in  interpolation  at  the  large  displacements. 

A convenient  failure  or  collapse  criterion  is  the  displacement 
beyond  which  the  structure  would  overturn  of  its  own  weight.  In  terms 
of  the  present  notation,  this  collapse  displacement,  Xc,  is  given  by. 
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= 1 + 1/0.014  = 72  for  the  3.4a  bridge  section. 

From  Fig.  5.16,  collapse  is  thus  assured  for2/ 

F = C.P,(0)  > 2.5R 
mad  p 

For  the  range  of  resistance  considered  in  the  response  computations  for 
Shot  9,  the  collapse  load  will  thus  vary  between  about  6 and  12  psi. 
Inasmuch  as  the  peak  load  for  Shot  10  is  estimated  to  exceed  27  psi 
(based  on  an  over-all  drag  coefficient  of  3)*  the  result  of  the  present 
analysis  is  entirely  consistent  with  the  actual  behavior  of  the  bridge. 
In  fact,  one  would  predict  failure  or  large  permanent  set  even  if  the 
over-all  drag  coefficient  were  taken  as  low  as  1. 


Table  5»1  - Results  of  Transient  Analysis 

at  at 

f(t)  = Fne  + F„e  + F.e  ^ (t  in  sec) 


Strain 

Record 


F2 

F3 

0.423 

-0.537 

Numerical  Constants 
of  Forcing  Function 


2 a3 

(l/sec)  (l/sec) 


-115 
-397 
• 8.51 


Frequencies 
of  Transfer 
Function 

w 1 
(cps) 

W2 

(cps) 

^5 

150 

48 

27 

37 

32 

120 

23 

91 

29 

71 

18 

62 

32 

81 

13 

66 

35 

39 

The  bridge  would  actually  be  expected  to  fail  at  some  lesser  dis- 
placement than  72  yield  displacements.  However,  in  the  region  of 
collpase  the  displacement  is  an  extremely  sensitive  function  of  peak 
load,  and  it  is  not  unrealistic  to  base  the  collapse  load  on  this 
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As  In  the  response  computations  for  Shot  9,  the  present  results 
can  be  given  no  further  meaning  other  than  to  state  that  the  simplified 
loading  and  response  schemes  described  above  provide  a realistic  esti- 
mate as  to  the  actual  behavior  of  the  3.4a  truss  bridge  in  Shot  10. 
Whether  this  agreement  is  mostly  coincidental  or  indeed  reflects  on  the 
general  validity  of  the  load  prediction  scheme  for  the  precursor  region, 
cannot  of  course  be  known  on  the  basis  of  one  result. 
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Fig.  3.1  Results  of  Transient  Analysis, Strain  Record  3-4aS2,  Shot  9 
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Fig.  5.2  Comparison  of  Experimental  Strain  Record  3.UaSl,  Shot  9,  with  Calculated  Strain  Record 
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Fig.  5.3  Results  of  Graphical  Analysis,  Strain  Record  3.4aSl,  Shot  9 


Fig.  5.4  Results  of  Graphical  Analysis,  Strain  Record  3.4aS2,  Shot  9 
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Fig.  5.8  Horizontal  Force  on  Truss  Bridge  Section  3.4a,  Shot  9 
(Based  on  Strain  Records  SI  and  S2) 
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Fig.  5.10  Assumed  Elastic-Plastic  Resistance  Function 
for  3«Ua  Truss  Bridge  Section 
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Fig.  5.H  Generalized  Resistance  Function 
for  Newmark-Brooks  Solution 
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Fig.  5.15  Approximate  Horizontal  Precursor  Loading, 
3.4a  Truss  Bridge  Section,  Shot  10 
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CHAPTER  6 


CONCLUSIONS  ARP  RECOMMENDATIONS 


6.1  CONCLUSIONS 

Most  of  the  strain  data  did  not  yield  Interpretable  results,  and 
therefore  none  of  the  stated  test  objectives  have  been  achieved.  In 
view  of  the  limited  and  uncertain  nature  of  the  experimental  data  ob- 
tained, it  is  very  doubtful  If  the  test  objectives  could  have  been 
realized  in  any  event.  But  the  fact  remains  that  it  was  not  possible 
to  interpret  adequately  the  data  which  were  obtained.  Thus,  no  posi- 
tive statement  can  be  made  even  with  regard  to  the  reliability  of  the 
data  Itself. 

Two  methods  of  analysis  were  attempted:  one  was  the  method  of 
transient  analysis,  and  the  other  was  essentially  a graphical  curve- 
fitting process.  A first  application  of  transient  analysis  proved  un- 
satisfactory due,  apparently,  to  an  inadequate  analytic  representation 
of  the  strain  data;  a second  application  gave  more  promising  results 
but  this  work  was  not  completed.  On  the  basis  of  present  knowledge  it 
is  not  possible  to  conclude  definitely  as  to  the  feasibility  of  transi- 
ent analysis  in  the  present  application. 

The  second  method  of  analysis,  which  consisted  of  graphically  sup- 
pressing the  obvious  harmonic  components  of  the  response,  provided  some 
qualitative  verification  of  the  pretest  load  prediction  method  for  the 
truss  bridge  section  but  did  not  yield  sufficient  quantitative  infor- 
mation to  satisfy  the  test  objectives.  As  a consequence  of  this  graphi- 
cal approach,  it  now  appears  that  a more  sophisticated  method  of  data 
reduction  along  the  lines  of  transient  analysis  is  mandatory  if  the 
desired  interpretation  is  to  be  made  of  the  response  of  a relatively 
complicated  dynamic  system.  In  this  connection  it  should  be  emphasized 
the  entire  subject  of  data  reduction,  as  it  pertains  to  the  problem  of 
net  force  measurement  systems,  deserves  a more  comprehensive  research 
effort  than  could  be  given  to  it  in  the  present  program. 

There  are  certain  lessons  to  be  learned  from  this  test  which  would 
apply  to  any  similar  data-reduction  scheme.  It  is  now  apparent,  and 
possibly  should  have  been  at  the  outset,  that  as  much  information  as 
possible  should  have  been  obtained  regarding  the  nature  of  the  dynamic 
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systems  being  tested.  That  is,  pre-  and  post- test  static  calibration 
tests  should  have  been  conducted  for  all  of  the  test  items.  This  would 
serve  to  determine  certain  elastic  properties  of  the  test  items  as  well 
as  to  provide  a check  on  the  strain  gage  system.  In  addition,  it  would 
have  been  most  desirable  to  conduct  suitable  dynamic  tests  in  order  to 
determine  the  principal  frequencies  and  damping  of  the  test  items.  In- 
formation of  this  type  would  have  undoubtedly  aided  in  the  study  and 
interpretation  of  the  strain  data. 

Perhaps  the  most  important  result  of  the  test  program  has  come 
about  through  consideration  of  the  damage  sustained  by  the  J.ka  truss 
bridge  section.  It  was  found  that  a simplified  dynamic  response  analy- 
sis, incorporating  the  (pretest)  predicted  loading  in  the  regular  re- 
flection region  for  Shot  9 and  a tentative  load  prediction  scheme  ap- 
plicable to  the  precursor  region  for  Shot  10  provided  an  adequate  esti- 
mate of  the  damage  sustained  by  the  bridge  in  both  instances.  While 
the  assumptions  of  the  response  analysis  are  suff iciently  uncertain  to 
invalidate  this  agreement  as  a check  on  the  essential  accuracy  of  the 
load  prediction  methods,  confidence  in  the  utilization  of  existency 
methods  for  damage  prediction  estimates  of  open-framed  structures  is 
certainly  increased  as  a result  of  this  test. 

6.2  RECOMMENDATIONS 

While  the  present  test  has  been  rather  unsuccessful,  it  is  believed 
that  the  original  objectives  are  still  valid  ones,  and  can  be  fully 
satisfied  only  by  means  of  full-scale  testing.  The  experience  of  this 
test  indicates  that  much  additional  work  needs  to  be  done  with  respect 
to  net  force  measurement  systems,  both  as  to  physical  design  and  inter- 
pretation of  data  before  such  tests  can  be  successfully  carried  out. 
Despite  the  present  unsuccessful  utilization  of  the  test  structures  as 
dynamometers  with  strain  output,  this  approach  has  obvious  advantages. 

It  is  recommended,  therefore,  that  a comprehensive  study  be  initiated 
on  the  data-reduction  problems  associated  with  net  force  measurement 
systems.  The  objectives  of  this  study  would  be  to  (a)  devise  methods 
of  analysis,  (b)  specify  the  class  of  physical  systems  whose  response 
will  successfully  permit  the  inverse  solution,  and  (c)  to  determine  the 
type  and  quality  of  the  necessary  data  to  permit  such  solutions.  It  is 
felt  that  particular  consideration  should  be  given  to  the  analog  computer 
technique  of  data  analysis  presented  briefly  in  this  report. 
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TRANSIENT  ANALYSIS 


A.l  INTRODUCTION 

In  the  present  application,  the  prime  data  (i.e.,  measured  strains) 
do  not  represent  in  a direct  sense  the  desired  information  (i.e., 
net  forces).  In  fact,  the  data  obtained  represent  the  response  or  a 
rather  complicated  dynamic  system  to  an  input  which  is  itself  the  de- 
sired information.  The  dynamic  system  consists  of  aerodynamic,  elastic, 
plastic,  electronic,  thermal  and  possibly  stochastic  components  coupled 
in  many  and  varied  fashions.  The  problem  at  hand  is  to  remove  from  the 
prime  data  the  influence  or  contributions  of  these  various  components 
and  of  the  system  itself,  leaving  the  desired  information.  Since  it  is 
difficult  (probably  impossible)  to  assess  the  exact  contribution  or 
influence  that  the  dynamic  system  contributes  to  the  prime  data,  a 
rational  approach  to  the  problem  is  to  attempt  to  assess  the  contribu- 
tion or  influence  of  those  components  which  could  be  expected  to  have 
the  most  pronounced  effects.  Due  to  lack  of  sufficient  knowledge  of  the 
test  conditions  and  the  reliability  and  characteristics  of  the  various 
recording  and  data-reduction  instruments,  the  magnitude  of  the  contri- 
butions of  some  of  the  above  components  cannot  even  be  estimated.  None- 
theless, it  has  been  assumed  that  the  major  influence  on  the  prime  data 
stems  from  an  elastic  response  of  the  cylinder-support  system. 

The  following  analysis,  referred  to  as  "transient  analysis,"  is 
designed  to  remove  the  influence  of  the  elastic  components  in  the  prime 
data,  which  is  undoubtedly  strongest  regardless  of  other  effects.  The 
analysis  as  applied  here,  represents  the  function  or  input  required  to 
produce  the  prime  data  as  the  response  of  an  elastic  model. 


A. 2 DERIVATION  OF  THE  BASIC  EQUATION 


be  effective  are  the  rigid-body  translation  and  rotation  of  the  founda- 
tion, rigid-body  motions  of  the  cylinder,  bending  and  axial  vibrations 
of  the  cylinder  supports,  bending  of  the  cylinder,  torsional  deflections 
of  the  cylinder,  or  vibrations  of  the  cylinder  as  a shell.  Regardless 
of  the  number  of  degrees  of  freedom  present,  the  equations  of  elastic 
motion  may  be  written  in  the  form: 
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fjtt) 


(A.l) 


j = 1,  2,  ...  N 

where  the  aij,  bj.j,  and  mj  are  real  constants,  and  f j(t)  is  the 
forcing  term  for  the  jth  degree  of  freedom.  X i(t)  is  the  response  of 
the  jth  degree  of  freedom  represented  by  the  jth  mass,  Equation 

A.l  is  the  most  general  system  of  equations  for  N degrees  of  freedom 


the  darap- 
which  may 


f <(t)  arise  pri- 


without  inertial  coupling.  The  constants  aji  and  b^j  are  the  damp- 
ing and  elastic  coupling  coefficients,  respectively,  some  of  which  may 
be  zero. 

For  the  present  application,  the  forcing  terms  f j(t)  arise  pri- 
marily from  the  pressure  acting  on  the  cylinder.  It  is  reasonable  to 
assume  that  motion  of  the  foundation,  if  it  occurs,  can  result  only  from 
coupling  with  the  cylinder  motion.  Thus,  only  those  of  Equation  A.l 
corresponding  to  masses  mj  associated  with  motion  of  the  cylinder  are 
non-homogeneou8  (i.e.,  fjft)  /0). 

The  masses  mj  associated  with  the  motion  of  the  cylinder  are  of 
two  types;  namely,  total  mass  corresponding  to  rigid-body  motions,  or 
generalized  modal  masses  for  the  bending  or  shell  vibrations.  The  forc- 
ing terms  for  each  of  these  result  from  either  the  integrated  pressure 
over  the  cylinder  alone  or  the  integrated  pressure  times  a modal  de- 
flection curve  over  the  cylinder.  If  the  pressure  is  assumed  to  be 
nearly  constant  in  the  axial  direction  along  the  cylinder,  these  forc- 
ing terms  will  be  essentially  proportional  to  one  another  at  each  in- 
stant of  time  so  that  their  time  histories  will  be  similar.  In  effect, 
Equation  A.l  will  contain  the  same  rlghthand  members  up  to  a constant 
factor,  which  may  be  zero. 

In  this  case.  Equation  A.l  becomes, 


”A(t)  + A aiA(t)  + , Vti(t) ' V(t) 


i =»  1 


i = 1 


(A.2) 


i = 1,  2,  ...  N 

It  is  well-known  (see  S.  Timoshenko,  Vibration  Problems  in  Engineering) 
that  the  response^  Xi(t) , Equation  A.l,  subject  to  homogeneous  initial 
conditions,  can  be*  given  by  a sum  of  Duhamel-type  integrals  as: 
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where  the  h^j  are  referred  to  as  transfer  functions  and,  for  a fixed 
(k, j) , represent  the  impulsive  response  in  the  jth  degree  of  freedom  to 
a unit  impulse  applied  in  the  kth  degree.  In  other  words,  hkj(t)  is 
the  solution,  Xj(t),  of  the  corresponding  homogeneous  system  of  equa- 
tions jj 

+ ^ b,  ^(t)  « 0 (A. 4) 

i = 1 


J = 1,  2,  ...  N 

subjected  to  the  inhomogeneous  initial  conditions 


Xj(0)  = Xj(0)  - 0 

Xk(0)  = 0,  X^O)  =1,  J + k 

Such  a solution  Xj(t)  * hk,(t)  is  known  to  be  a linear  combination  of 
exponential  functions  ( in  general , complex)  A/  of  the  form  ( see  S . Timo- 
shenko, Vibration  Problems  in  Engineering) . 

^ a“J'Xnt  <A,5) 

where  * 0,  + 1,  ...  + M ) are  the  2B  roots  of  the  characteristic 

equation  (system  of  equations) 

N N 

"•> x2  + i ?i  *1J  X +i?iblJ"° 

Jsl,  2,  • • • N 

If  Equation  A. 5 is  substituted  in  Equation  A. 3,  there  results: 

“ r'  ” k1  xft-n 

X-,(t)  « 2.  J tA  T)  £ o*Je  4T(a.6) 

j k = 1 o K n = -N  n 

If,  in  addition,  the  fk(t)  are  timewise  proportional  so  that 

ffc(t)  - ckf(t)  (A. 6a) 

as  in  Equation  A. 2, then  Equation  A. 6 may  be  written, 

1/  It  is  assumed  here  that  the  2N  characteristic  roots  are  distinct. 
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J o k=ln=-N*n 


Xj(t) 


f(T)  2 ^ e“" 


n - -N 


XJt  -T  ) 


dr 


XJ(t) 


f(r)  hj(t  - t)  dT 


(A. 7) 


where: 

i 

n 


It  is  seen  that  the  quantity  hj  (t  - r ) in  Equation  A. 7 is  precisely 
the  form  of  the  transfer  functions  h^jCt)  in  Equation  A. 3 (independent 
of  k,  however)  and  is  a linear  combination  of  exponential  functions  in- 
volving only  the  characteristic  roots  of  the  system.  Equation  A. 7 is 
the  desired  basic  equation  and  may  be  used  to  obtain  the  unknown  forc- 
ing function  f(t)  from  a known  response  Xj(t).  The  quantity  f(t) 
actually  represents  only  the  shape  of  the  forcing  function,  since  it  is 
defined  up  to  the  constant  factors  c j of  Equation  A. 2. 

The  problem  of  uniqueness  and  the  effect  of  "errors"  in  the  response 
curve  Xj(t)  &s  they  affect  the  solution  for  f(t)  obtained,  will  be 
discussed  at  the  appropriate  steps  in  the  subsequent  analysis.  It  will 
be  seen  that  the  uniqueness  problem  essentially  reduces  to  a problem  in 
curve-fitting  and  may  be  treated  in  several  ways. 


A. 3 SOLUTION  OF  THE  BASIC  EQUATION 


Equation  A. 7 can  be  solved  by  a variety  of  methods,  none  of  which, 
however,  leads  to  a unique  solution  for  the  forcing  function  f(t).  For 
the  present  application  it  is  assumed  that  f(t)  can  be  represented  by 
exponential  functions  in  the  form: 
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f(t)  = FLe 
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(A. 8) 


where  the  and  ai  are  constants.  The  number  (K)  of  terms  in 

Equation  A. 8 will  be  determined  by  the  complexity  of  the  actual  forcing 
function  or,  at  least,  the  complexity  required  of  f(t)  in  order  to 
produce  the  known  response  X(t).  It  should  be  mentioned  that  the  form 
of  Equation  A. 8 in  no  way  restricts  the  generality  of  application. 

If  Equation  A. 8 is  substituted  into  Equation  A. 7,  and  Equation 
A. 7a  is  used,  there  results  upon  integration: 
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where, 


ai  = 
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n = -N 
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(A.10) 
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It  is  assumed  here  that  a±  ft  A 


J 

for  each 


i and  n,  since  other- 


wise a resonance  condition  would  result,  which  certainly  does  not  happen 


in  the  physical  case  considered  here.  Equation  A. 9 gives  X*(t)  as  a 


sum  of  exponential  functions,  and  the  present  problem  reduces  to  that 
of  fitting  such  a sum  to  the  prime  data  representing  Xj(t)  for  a dis- 
crete set  of  times,  t. 

A curve-fitting  scheme  results  in  the  constants,  A^,  Gt^,  Bn,  and 
An  (real  and  complex,  in  general)  which,  in  turn,  may  be  used  to  com- 
pute the  desired  Fj[  through  Equation  A. 10.  Actually,  Equation  A. 10 
wilx  also  yield  the  fi  “ which,  however,  are  not  of  direct  concern  in 


the  present  application.  With  the  Gj_,  A^,  Bn>  and  Xn  known,  Equa- 


tion A. 10  constitutes  a system  of  2N  + K equations  in  the  2N  + K 


unknowns, fi  £ and  Fi.  With  the  and  F<  determined,  the  desired 

forcing  function  is  then  given  by  Equation  A. 8.  Actually,  a little  more 


information  is  assumed  here;  namely,  that  not  only  the  constants  Ct 


Ai>  Bn, 


and  X 


of 


n are  obtained  from  the  curve-fitting,  but  that  the  sets 
Ai's  and  aq's  can  be  distinguished  from  the  set  of  Bn's  and 


1 


n's.  Otherwise,  the  structure  of  Equation  A. 10  would  not  be  determined 


so  that  a unique  forcing  function  could  not  be  constructed  through  Equa- 
tion A. 8.  In  the  present  application,  the  distinction  between  the  two 
classes  of  terms  in  Equation  A. 9 can  be  made  through  the  exponents  a i 
and  X n themselves.  Since  Equation  A.l  can  be  assumed  to  be  stable. 


none  of  the  characteristic  roots 


can  be  pure  real.  On  the  other 
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hand,  it  is  assumed  that  the  forcing  function,  f(t),  contains  no  har- 
monic components  so  that  a realistic  curve-fitting  scheme  should  yield 
only  real  exponentials  for  f(t)  in  Equation  A. 8.  In  brief,  the 

non- real  exponents  obtained  from  curve-fitting  are  associated  with  the 
transfer  function 


Mt) 


N 

1 


-N 


and  give  the  a n,  while  any  real  exponents  obtained  are  associated  with 
the  forcing  function 


f(t) 


The  appropriateness  of  this  choice  of  association  can  be  asserted  only 
after  some  study  is  made  of  the  effect  of  errors , etc . , of  the  prime 
data  on  the  computed  exponentials.  It  would  be  desirable,  of  course, 
to  obtain  the  natural  frequencies  of  the  system  either  analytically  or 
experimentally  independent  of  the  above  analysis,  since  then  they  could 
be  identified  in  Equation  A. 9, and  f(t)  could  then  contain  harmonics 
without  affecting  the  above  solution. 


A.k  CURVE-FITTING 


Equation  A. 9 may  be  written  in  the  form: 
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X(t)  = £ 


k = 1 


Cke 
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(A. 11) 


where  the  and  7 ^ are  in  general  complex  numbers,  and  where  the  sub- 
script j has  been  dropped,  since  the  only  concern  is  with  the  response 
of  but  one  component  of  the  system.  Curve-fitting  consists  of  finding 
those  values  for  c^  and  7 k such  that  X(t)  as  given  by  Equation  A. 11 
best  fits  the  given  prime  data. 

Method  A.  The  most  direct  method  of  curve-fitting  is  accomplished 
by  visual  inspection  of  plots  of  the  prime  data  itself.  A discussion 
of  this  method  and  the  results  of  its  use  in  the  present  application  are 
presented  in  Section  5.1.2. 

Method  B.  A second  method  consists  of  rather  standard  technique 
from  the  theory  contained  in  the  calculus  of  finite  differences  (see 
The  Calculus  of  Observations,  Whittaker,  E. , and  Robinson,  A.).  It  is 
assumed  now  that  the  prime  data  are  given  in  equal  interval  tabular  form. 

Let  the  interval  of  tabulation  of  the  function  X(t)  be  A , and 
designate  X(t)  at  the  time  when 


t = n A 


n = 0,  1,  . . . P 
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as  Xn. 
as 


The  assumed  form  of  X(t) (Equation  A. 11)  may  then  be  written 


(A. 12) 


where 


(A. 12a) 


This  method  of  curve-fitting  rests  upon  the  observation  that  Equation 
A. 12  is  the  solution  to  a linear  finite  difference  equation  with  con- 
stant coefficients.  This  finite  difference  equation  may  be  written  for 
each  consecutive  group  of  (L  + 1)  tabulated  values  of  Xn  as 

L 

E D/  Xn  / = 0 (A.13) 

X = o * n+  1 

where  = 1 and  the  remaining  D's  are  to  be  determined.  If  (P  + 1) 
is  the  total  number  of  tabulated  values  of  Xn,  Equation  A. 13  may  be 
written  (P  - L + l)  times,  each  equation  corresponding  tc  each  of  the 
values  of  n in  the  range 


n - 0,  1,  2,  ...  (P  - L) 

Consider  the  resulting  set  of  simultaneous  equations 


DoXo  + D1X1  +,,,+DL-lXL-l-'XL 


_XL  + 1 


where  the  X^s  are  known  from  the  tabulated  function  X(t).  This  set 
of  simultaneous  equations  consists  of  (P  - L + 1)  equations  in  the 
(L  * l)  unknown  D's.  A solution  to  the  system  of  equations  is  possible 
if  P > 2L.  When  L,  the  number  of  exponential  terms  used  to  approxi- 
mate  X(t),  is  chosen  equal  to  P/2,  the  system  has  a unique  solution 
if  the  determinant  of  the  coefficients  of  Xn  is  non-zero.  Normally, 
this  condition  exists  and  the  unknown  D's  are  determined  by  ordinary 
methods  (see  The  Calculus  of  Observations,  by  Whittaker  E.,  and  Robin- 
son, A.)  and  are  used  to  form  the  characteristic  equation  which  is  the 
polynomial 

0 = Do+D^X+DgX+»«»  + Dk_^X  +X  (A.  15) 


(A.14) 
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The  roots  of  Equation  A.15  are  the  X k's  referred  to  in  Equation  A. 12 
and,  hence,  may  be  used  to  determine  the  /it'3  which,  in  turn,  are  used 
to  form  the  approximate  representation  of  X(t). 

The  actual  determination  of  the  yk's  is  considered  now,  since 
the  roots  of  Equation  A.15  may  be  positive,  negative,  or  complex.  Re- 
gardless of  the  type  of  root,  X k can  be  written  as 


X 


k 


(i  = 


) 


where  9 is  equal  to  zero  or  pi  ( T ) if 
respectively.  If  X ^ is  a complex  root, 
angle  and  lies  between  - if  and  7T  . if 


xk  is  positive  or  negative, 
9 is  the  Cauchy  principal 
Xk  is  positive, 


if  xk 


log 


is  negative. 


(A. l6a) 


rk = a"  108  ('V  + X (A,l6b) 

In  the  latter  case  the  imaginary  part  of  Xk  is  dropped,  since  it  in- 
troduces an  extraneous  frequency  whose  period  is  2 A . The  conclusion 
that  this  frequency  is  extraneous  is  based  upon  the  observation  that  it 
is  dependent  only  upon  the  choice  of  the  tabulation  interval,  and, 
therefore  may  be  changed  to  any  value  by  a suitable  choice  of  interval 
Thus,  for  all  real  values  of  Xk 


= IT  108 


(A.17a) 


£7  An  alternate  method  of  dropping  the  frequency  introduced  by  the  nega- 
tive roots  is  to  consider  only  the  even  terms  in  Equation  A. 12.  Equa- 
tion A. 12  becomes 


2n 


n 


And  Equation  A.  12a  b. 


from  which 


>e  comes 

%(2A) 


r*- 


aS~  108  (Xk’ 


and  X(t)  again  has  the  form  of  Equation  A. 12  where  the  extraneous  fre- 
quency corresponding  to  W /A  is  not  present. 
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If  is  complex,  it  must  be  one  of  a conjugate  pair  and  both  pairs 

are  contained  in 


(A. 17b) 


Having  obtained  the  Xk's  from  Equations  A. 17a  and  A. 17b,  the  c^'s 
may  be  obtained  from  Equation  A. 12  in  combination  with  the  initial  con- 
ditions. These  conditions  are  satisfied  if 


M 

2 c = 0 (A.18) 

k = 1 K 


and 
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k = 1 


\ =*  - 0 


(A.19) 


Since  Equation  A. 12  is  to  be  satisfied  for  each  tabulated  value  of  X(t), 
there  results  the  following  P equations  in  (L  - 2)  unknowns, 


rlXl+  C2X2+  **•  + CLXL  = X1 

2 2 2 
C1  X1  + c2  X2  + **•  + CLXL  = X2 

C1X1+C2X2+  ”•  CLXL  = XP 


(A. 20) 


where  Equations  A.18  and  A. 19  are  to  be  used  to  eliminate  two  of  the 
cjt's.  The  above  system  of  equations  is  overdeterminate  and  usually  in- 
consistent, since  a finite  number  of  terms  is  used  to  approximate  X(t). 
The  method  of  least  squares  is  used  to  determine  the  most  probable  set 
of  (L  - 2)  equations  in  the  (L  - 2)  unknowns,  (see  The  Calculus  of 
Observations . Whittaker ,E . , and  Robinson,  A.)  and  then  this  set  of  equa- 
tions is  solved  for  the  c^'s.^/ 

Method  C.  A third  method  of  curve-fitting  utilizes  an  analog  com- 
puter, and  represents  a rapid  trial-and-error  process.  A discussion  of 
this  method  is  presented  in  Section  5.1.^. 


A. 5 DISCUSSION  OF  APPLICATION  OF  TRANSIENT  ANALYSIS  (Method  B) 


The  application  of  transient  analysis  as  described  in  Method  B 
consists  of  the  following  steps: 

a.  The  frequencies  and  exponents  occurring  in  the  expansion  of  the 
response  curves  (prime  data)  are  determined  by  analogy  with 

J?  If  L < P/2,  the  method  of  least  squares  is  used  to  reduce  the  equa- 
tions to  a determinate  number . 
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finite  difference  equations.  The  prime  data  are  interpreted 
as  representing  the  solution  of  a system  of  linear  finite 
difference  equations.  The  characteristic  equation  (polynomial) 
for  this  system  of  difference  equations  is  constructed  using 
the  prime  data  and  solved  to  obtain  the  characteristic  roots. 
These  roots,  in  turn,  are  taken  as  the  desired  frequencies  and 
exponents  of  the  response  curve  expansion. 

b.  Having  obtained  the  frequencies  and  exponents  of  the  response 
curve  expansion,  the  amplitudes  associated  with  each  of  the 
exponential  terms  are  then  determined  by  a straightforward 
least  squares  fit  to  the  prime  data. 


c.  With  the  expansion  of  the  response  curve  determined,  the  forc- 
ing function  is  then  associated  with  the  real  exponential  terms 
and  is  obtained  through  Equation  A. 10. 

The  success  of  the  application  of  transient  analysis  to  the  pre- 
dicting of  actual  forcing  functions  hinges  upon  several  reasonable  but 
crucial  assumptions.  First,  and  of  utmost  importance,  is  the  assump- 
tion that  the  dynamic  system  producing  the  response  curves  be  essenti- 
ally linear.  Just  what  magnitude  of  non-linear  effects  such  as  stochas- 
tic errors,  or  baseline  shifts  are  tolerable  is  not  known  at  present, 
and,  indeed,  represents  a major  area  of  required  investigation  before 
any  net  force  measurement  technique  of  the  present  type  can  be  evaluated. 
The  effect  of  errors  (known  or  unknown)  in  the  prime  data  affect  the 
analysis  most  critically  in  step  a above.  In  fact,  if  the  frequencies 
and  exponents  determined  in  step  a are  realistic,  the  success  of  the 
entire  analysis  appears  to  rest  only  upon  the  identification  made  in 
step  c,  which  will  be  discusseu  presently.  One  criticism  of  the  pro- 
cedure in  step  a for  determining  the  all-important  frequencies  and  ex- 
ponents appears  immediately  as  the  result  of  using  the  prime  data  it- 
self. Since  these  data  in  the  present  application  are  known  to  contain 
various  types  of  errors,  it  would  be  desirable  to  incorporate  some  data- 
smoothing  process  prior  to  the  application  of  transient  analysis.  There 
are  several  types  of  data- smoothing  techniques  available  with  various 
degrees  of  reliability  and  meaningfulness  to  the  present  application, 
but  the  most  promising  appears  to  be  the  use  of  the  autocorrelation 
function.  For  the  latter  application  step  a would  use  the  autocorrela- 
tion function  defined  as 


The  autocorrelation  function  X(t)  contains  the  same  frequencies 
and  exponentials  as  do  the  prime  data  X(t)  and  smooths  the  data  by 
eliminating  stochastic  and/or  bias  errors.  This  smoothing  technique  is 
being  tried  and  tested  at  present  in  other  and  similar  applications 
(see  Operation  TEAPOT  Project  3.2,  Study  of  Drag  Loading  on  Structures 
in  and  out  of  the  Precursor  Zone,  WT-1124) . 

Of  course,  the  analogy  technique  applied  in  step  a does  itself 
perform  some  sort  of  averaging  on  the  data  fed  in.  Just  what  the  nature 
of  this  averaging  is,  is  not  known  at  present.  A perturbation  study  of 
the  effect  on  the  exponents  and  frequencies  as  determined  by  step  a is 
necessary  before  a realistic  evaluation  of  transient  analysis  can  be 
made.  Because  of  the  complications  involved,  such  a study  would,  no 
doubt,  have  to  be  carried  on  in  the  sphere  of  numerical  experimentation. 

Step  b determines  a "best  fit"  of  the  prime  data.  The  fit  obtained 
itself  is,  in  many  respects,  a measure  of  the  reliability  of  transient 
analysis  in  the  present  application.  However,  a realistic  criteria  for 
"goodness"  of  fit  has  not  yet  been  established;  at  least,  not  in  terms 
of  the  reliability  or  meaningfulness  of  the  forcing  function  so  obtained. 

A point-by-point  check  on  the  fit  may  not  be  too  meaningful,  since  the 
analysis,  to  be  successful,  must,  of  necessity,  ignore  "errors"  in  the 
prime  data  and  pick  out  only  pertinent  information.  However,  the  method 
of  least  squares  as  applied  in  step  b appears,  at  present,  to  have  the 
most  reliable  characteristics  for  curve-fitting  as  applied  in  this  ap- 
plication. 

The  procedure  of  step  c is  based  on  the  assumption  that  the  forc- 
ing function  is  represented,  or  "sufficiently  well"  represented,  by  a sum 
of  real  exponentials.  Just  what  is  "sufficiently  well"  to  assure  suc- 
cess of  transient  analysis  has  not  yet  been  determined.  The  limits  of 
reliability  could  most  likely  (and  probably  of  necessity)  be  determined 
by  numerical  experimentation,  also. 

Finally,  the  assumption  that  all  forcing  functions  are  timewise 
proportional  needs  to  be  considered.  The  first  question  to  be  answered 
is  whether  the  forcing  functions  considered  are  essentially  proportional 
or  not.  The  second  question  to  be  answered  is  what  is  "sufficiently 
proportional"  in  terms  of  success  of  the  application  of  transient  analysis. 
Thirdly  (since  the  first  two  questions  have  not  as  yet  been  answered), 
what  happens  to  transient  analysis  if  the  forcing  functions  are  not  es- 
sentially timewise  proportional.  In  particular,  does  the  forcing  func- 
tion obtained  by  applying  transient  analysis  in  this  latter  case  have 
any  significance  or  yield  any  real  information.  Of  these  three  ques- 
tions, only  the  third  can  be  answered  (to  some  degree  at  least)  at  the 
present  time. 

In  answer  to  the  question  of  what  happens  if  the  forcing  functions 
are  not  timewise  proportional,  the  analysis  of  the  preceding  section, 
is  repeated  in  brief  as  follows: 

Equation  A. 6 is  written  in  the  form 

N -t  N 

X(t)  = 'Z  J f.(T)  2 
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(A. 22) 


UNCLASSIFIED 


86 


where  the  index  j has  been  dropped,  since  the  only  concern  is  the 
response  in  a single  degree  of  freedom.  The  forcing  function  fj^t) 
is  written  in  the  form 

I v a.t 

fk(t)  = 2 F k e 1 (A. 23) 

i = 1 1 

where  the  ai  range  over  all  exponentials  appearing  in  any  of  the 
degrees  of  freedom,  so  that  some  of  the  may  be  zero.  With  Equa- 

tion A. 23  substituted  in  Equation  A. 22  and  the  integration  performed, 
there  results  the  equation: 
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(A. 24) 


i = 1 * n * -N 

which  is  precisely  the  same  form  as  Equation  A.Q.  However,  the  and 
Bn  are  related  to  the  forcing  coefficients  F^*  and  the  transfer  func- 
tion coefficients  a n*  by 
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Bn  = k?i  i ^l  ai  ‘ Xn 


which  are  the  counterparts  of  Equation  A. 10.  The  relationships  between 
Equation  A. 10  and  A.25  results  in  (dropping  the  index  j) 


N 


2 F k a k = $ F 

k-1  1 n Dl 


(A. 26) 


n = + 1,  + 2,  . . . + N 
i = 1,  2,  ...  K 

Hence,  the  curve-fitting  of  step  b may  be  applied  to  Equation  A. 24  re- 
sulting in  the  aj_,  a B~  and  X n which,  in  turn,  are  used  in  Equa- 

tion A. 10  to  determine  the  p n and  Fj_.  However,  step  c is  not,  in 
general,  valid  since  the  Fj_  are  not  the  actual  forcing  function  coef- 
ficients. In  fact,  there  is  no  longer  a single  forcing  function  in 
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question,  but  one  for  each  degree  of  freedom.  The  actual  forcing  func- 
tion coefficients  Fik  are  related  to  the  /3n  and  F^  by  Equation 
A. 2 6.  However,  in  Equation  A. 26  only  the  J3  n and  Fj^  are  known,  while 
both  sets  of  Fj_k  and  a nk  are  unknown.  Since  there  are  (2N)  x K 
equations  in  Equation  A. 26  and  (2N)  x K x N unknowns.  Equation  A. 26  doe6 
not,  in  general,  determine  the  actual  forcing  function  coefficients  F^ 
uniquely  and  an  infinitude  of  solutions  may  be  constructed.  It  should  be 
mentioned,  however,  that  some  information  about  the  forcing  functions 
has  already  been  determined  in  step  b.;  namely,  all  the  exponents  a ji 
appearing  in  any  of  forcing  functions  f^(t)  have  been  determined,  and 
Equation  A. 26  represents  a number  (even  though,  in  general,  not  enough) 
of  necessary  conditions  on  the  coefficients  F^.  Certain  auxiliary 
information  about  either  or  both  the  dynamic  system  or  forcing  functions 
obtained  analytically  or  experimentally  together  with  Equation  A. 26  can, 
in  certain  instances,  determine  a unique  set  of  coefficients  Fj_^.  For 
example,  if  curve-fitting  results  in  but  one  real  exponential,  transient 
analysis  yields  a unique  forcing  function  (in  a trial  way,  of  course). 

It  is  interesting  to  see  how  the  assumption  of  timewise  proportion- 
ality of  forcing  functions,  together  with  Equation  A. 26  yields  a unique 
solution  (as  it  must,  by  previous  considerations). 

The  assumption  of  proportionality  is  equivalent  to  Equations  A. 24, 

A. 6a,  and  A. 8: 


F . = c.  F . for  each  i and  k 

i k i 

so  that  Equation  A. 26  may  be  written 
N f N 
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(A. 28) 


n = + 1,  + 2,  . . . + N 
i 3 1,  2,  ...  K 

which  identifies  the  F±,  as  well  as  the 

£n  = £ ck  a l 

k = 1 

of  Equation  A . 'JeSJ  uniquely,  since  the  right  members  of  Equation  A. 28 
are  known.  In  effect  the  assumption  of  proportionality  uncouples  Equa- 
tion A. 26  (which  is  generally  coupled  in  the  independent  indices  i 
and  n)  into  two  independent  sets  of  equations 


Fi  " Fi  , 


i - 1,  2,  ...  K 


Dropping  the  index  j. 


UNCLASSIFIED 


and 


N 


1 

k = 1 


n = + 1,  + 2,  . • . + N 


Less  drastic  assumptions  regarding  the  forcing  functions  or  transfer 
functions  could  be  expected  to  yield  additional  information  about  the 
forcing  functions.  For  example,  if  Equation  A. 26  is  summed  with  respect 
to  i,  there  results 


N 

r k 1 

r k i 

2 

k = 1 

i = 1 

a k = /3 

n ^n 

2 

i = 1 1 

J 

(A.29) 


The 


inner  sum 

/ K 


2 

i = 1 


is  clearly  (see  Equation  A. 23)  the  quantity  ffc(O)  which  is  the  mag- 
nitude of  the  forcing  function  in  the  kth  component  at  time  t = 0. 
Hence,  Equation  A.29  may  be  written: 


N 


fk(0) 


(*.30) 


n = + 1,  + 2,  ...  + N 

It  is  possible  that  at  least  the  relative  magnitudes  (if  not  the 
absolute  magnitudes)  of  the  ffc(0)  can  be  estimated  from  theory,  etc., 
so  that  Equation  A.30  could  be  used  to  yield  additional  independent 
conditions  on  the  unknowns  of  Equation  A. 26.  If  the  absolute  magnitudes 
of  the  Fj^(0)  are  known.  Equation  A.30  yields  2N  additional  equations. 
If  only  their  relative  magnitudes  are  known,  2N-1  additional  equations 
result. 

Another  possibility  is  that  some  of  the  F-j.k  may  be  known  from 
other  considerations,  in  effect,  reducing  the  number  of  unknowns  in 
Equation  A. 26.  For  example,  some  degrees  of  freedom  may  not  admit  a 
forcing  function  and  are  known  to  be  excited  only  through  coupling  with 
other  degrees.  Hence,  the  corresponding  Fj_k  are  zero.  Certain  time- 
wise  proportionality  conditions  may  also  exist  between  certain  degrees 
of  freedom,  such  as  the  set  of  all  bending  modes,  further  reducing  the 
number  of  unknowns. 
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